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ABSTRACT 


In data communication systems periodic data streams occur 
due to (a) transmission of periodic sequences during idle condi- 
tion for ensuring efficient timing recovery , (b) transmission of 
training sequences where adaptive equalization is incorporated at 
the receiver and (c) periodicity in the message itself. The 
presence of periodic sequences may lead to certain difficulties 
which are dependant either on the period of the data or on 
pattern thereof. These are. (a) single tone interference in 
adjacent channel produced because of the periodic data streams, 

(b) inefficient adaptive equalization resulting due to short 
periodic data streams and (c) inefficient timing recovery resul- 
ting either due to infrequent changes in amplitude of the data or 
due to absence of densely spaced line spectrum during the idle 
period, depending upon the type of timing recovery scheme adapted. 
A device called "scrambler" is used to circumvent them. 

In this thesis we study the performance of binary digital 
data scramblers for periodic input sequences with respect to the 
period of the output sequence and spectral coefficients of input 
and output sequences. The periodic analysis of the output of the 
scrambler is based on D-transform techniques using the properties 
of polynomials and partial fraction expansion of partial fraction 
expansion of rational functions of polynomials over GF(2.). The 
spectral analysis of input and output sequences is based on 
spectral theory of linear sequential circuits over finite fields. 
It is shown that (a) if the period s of the input sequence is 
not a multiple of the period q of zero input response of the 



scrambler, the output period is of period equal to the least 

common multiple of s and q except for a unique state for which 

the output period is s itself and (b) if the period s is a 

multiple of q then the period of the output sequence is also a 

multiple of q. It is observed that for a n-stage scrambler 

(a) eigenvalues for all eigensequences over an extension field 
nr 

GF{2 “) is defined if n is not a divisor of m; otherwise eigen- 
values are not defined for n eigensequences, (b) spectral 
coefficients corresponding to eigensequences for which eigen- 
values are not defined, are zero for those input sequences of 
period s (s equal to or a divisor of 2 n -l); the new spectral 
coefficients, if any, appearing in the output sequence at 
these locations are due to the zero input response for a state 
determined by, the unique states corresponding to eigensequ- 
ences having non-zero coefficients at the input and initial 
state of the scrambler, and (c) for other input and output 
sequences both having periodicities equal to or a divisor of 
2 X -1 , i is not equal to n, no new spectral coefficients appear 
at output other than those determined by input coefficients 
and respective eigenvalues. 



CHAPTER 1 


INTRODUCTION 

In data communication systems periodic data streams 
may occur due to 

(a) transmission of periodic sequences during idle condition 
for ensuring efficient timing recovery, 

(b) transmission of training sequences where adaptive equali- 
zation is incorporated at the receiver, and 

(c) periodicity in the message itself. 

The presence of periodic sequences may lead to certain diffi- 
culties which are dependant either on the period of the data 
or on pattern thereof. These are, 

(a) Single tone interference- in adjacent channel produced * 
because of the periodic data streams. 

(b) Inefficient adaptive equalization resulting due to short 
periodic data streams. 

(c) Inefficient timing recovery resulting either due to 
infrequent changes in amplitude of the data or due to 
absence of densely spaced line spectrum during the idle 
period, depending upon the type of timing recovery scheme 
adapted . 

In this chapter we illustrate these difficulties and 
describe a device called "Scrambler" which is used to circum- 


vent them 
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1.1 Spectrum of Periodic Data 

Consider a gating function as shown in Figure 1.1. 

For the ease of illustration let us assume that the period is 
some multiple of the pulse width . Expanding this function 
using exponential Fourier series , its spectral components are 
given by 

F n = r sa <s r> <ia) 

where n = ..., -2, -1, 0, 1, 2, ... and Sa(x) is sampling 
function. The frequency spectrum is shewn in Figure 1.2. The 
fundamental frequency is 2lyT . The frequency spectrum is a 
discrete function and exists only at “ = 0, + ^- , + * ••• 

and has amplitudes — , ~ Sa(;~) / — Sa(— —) / ... respectively. 

Now consider a periodic data . . . , 1001 , 1001 , . . . as 
shown in Figure 1.3. This can be considered as made up of a 
gating function f(t) and a shifted gating function f(t - 36). 
The spectral components of f(t - 3 6) are given by 


F n 


AS 

T 


Sa 




-j a) (36) 


^ sa (2|i) 



( 1 . 2 ) 


The discrete frequency components which appear in f (t - 36) 
are same as those in f(t). If we denote the spectral components 
of the sequence . . . , 1001 , 1001 , . . . by F^ then 


F n 


n 


+ F' 


n 


(1.3) 


Amplitudes of any discrete frequency may get added up or get 
reduced since the phase angle is involved. There is atleast 
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one frequency component i.e. DC component the amplitude of which 
always gets added up (of course this is not the case in bipolar 
signalling) . In any case, amplitudes are inversely propor- 
tional to period T. Also note that as T increases/ the spectrum 
becomes densely spaced. 

1.2 Single Tone Interference 

New consider a typical basic group modulation scheme 
in a carrier communication system shown in Figure 1.4 [16], 

Let us assume that channel 1 has a device which introduces 
square law nonlinearity before modulation taking place. A tone 
of 2.5 KHz (say) in this channel gives rise to its second 
harmonic i.e. a 5 KHz tone. This 5 KHz tone after modulation 
appears at 103 KHz, i.e. in the domain of channel 2 (100 KHz - 
104 KHz) and causes interference. This is called as 'Single 
Tone Interference'. For this reason, system engineers place 
limits on the levels of isolated tones in a customer's trans- 
mission band. From previous section it is clear that tones 
are produced in data transmission system by periodic data 
sequences. Because of the inverse relation, the upper bound 
on tone levels is then gets translated to a lower bound on the 
period of the data sequence. Hence it is desirable to map any 
source sequence on to a periodic channel sequence of larger 
period. 

1.3 Adaptive Equalization 

Consider the block diagram of a communication system 
shown in Figure 1.5. For distortionless transmission of signal, 
it must experience a specified constant amount of gain or loss 
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Transimltter 



Figure 1.5 Block Diagram of a Communication System. 


and a linear phase shift i.e. 

jut 

A 


S(w) 


= K s(») e 


(1.4) 


where K and t are constants. This will be satisfied only if 
o 


X (») = T(«) H(«) R(w) = K e 


jut. 


(1.5) 


But in real systems, this is not the case as X e< ^(u) is given by 


X eq(“) = K ( w > e 


j0(u) 


( 1 . 6 ) 


This deviation from the ideal transmission characteristic 
causes distortion of input signal which causes overlapping 
between successive symbols in data communication and is known 
as 'Intersymbol Interference' (ISI). This can be combated by 
using a device called equalizer as indicated in Figure 1.6. 

Equalizer 



Jut, 


Figure 1.6. Block Diagram of a Communication System with 
Equalizer. 
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The transfer function E(o>) of the equalizer should be 

-j 0(“) - »(t -t*) 

E(») = T T -~r— Y e ° , for all w (1.7) 

W ' 

z\ Jdt 

Then S(w) = S(») X (to) E(») = K S(w) e (1.8) 

eq 


At high data rates — 4800 bits/sec and above — where 
the channel characteristics vary slowly compared to data rate , 
the setting of the equalizer is automatically adjusted to match 
the channel. This is known as 'adaptive equalization'. In 
the next section , we will see how the equalization gets degraded 
because of the presence of short periodic sequences in the data. 

1 .4 Spectrum of the Received Signal [5 3 

Now consider a digital transmission scheme illustrated 
in the Figure 1.7. Assuming an ideal channel, the signal 



Figure' 1.7 Block Diagram of a Digital Transmission Scheme. 
s(t) received at the receiver is given by 


s — t 

s(t) = > a g (t - nT) cos « t 

^ n c 

n 


(1.9) 


Let l a n> ^ P^i 0 ^ 0 with the period N. 
given by 


The DFT Of 



-jkn QT 

A (k S2 ) = y a e k = 0,1, ... , N“1 

n=0' 


( 1 . 10 ) 
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and the inverse relation is 


22 A(k 0 ) e jknn T , n = 0,1/.. , N-l 
k=0 


( 1 . 11 ) 


1 2n 1 

where T is the symbol duration and 0 = (— ) (~) = ^ “ 


N s 


where oi is symbol frequency in rad/sec. Hence the DFT has N 
s 

1 

components uniformly spaced ^ Hz apart and the spectrum repeats 

• i 

every — Hz. The spectrum S(&) of s(t) is given by ; 

S(w) = % a e “- ) " a5nT G(to) * 6(u-u .) , w > 0 (1.12) 

V r\ p 


and using Equation (1.11) in Equation (1.12) we obtain. 


S( ) 


i A(k Q ) e jkn Q T e~i “ nT G(u> ) * 6 ( w- “ c > 


= n . A(k } 5Z 

k=0 in 


-j(<a-ka)nT 


G(u) *6(“-w^_) 


(1.13) 


Using the identity / , e 

n 

(1.13) gives. 


-jn»T _ (la, _ 2 jEj in nation 

n 4 


S(w) = 


A(k& ) 


6( w - kfl - na> ) G(u) * 6( a - u ) .. 


i.e. 


s(») 


„ 2 NrL r 

> ] A(kfi ) > G (k n -f naj g ) 5( u - a) c - k n - n“ g ) 


.2 N-: 


S(M) = 



A(k n ) > ’G(k n+ nu) 6( B - to - kfl - nw ) 

s S C o 

n 


(1.14) 
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It is clear that S(co) has discrete tones at u c + kn+ n^. In 

particular data transmission systems with pulses of less than 

100% excess band width, G(k2+ nto ) will be zero if n is not 

s 

equal to zero or -1. Hence for n = 0 and -1, 



Figure 1.8 Spectrum given in Equation (1.15). 

The spectrum is illustrated in Figure 1.8. It is clear that 
the envelope A(k 8 ) modulates the base band pulse shape 
G(w - a) ), in the range w - to to to + to . 

C C o G S 

In adaptive equalization, since s(t) is used to adjust 
the tap weights, it is obvious that a constant A(k ft ) envelope 
will be ideal. A(k ft ) has frequency components at discrete 
points and hence the equalization can be achieved only at those 
discrete frequencies which are in the range of - w s to a> s . For 
example if the period were two symbol intervals (i .e . 2T ) , from 
Figure 1,2 it is clear that there are only two discrete frequency 
components one at - -sp and the other at + • in the Nyquist 
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band. Hence the equalizer can correct for distortion only at 
these two points. Consequently at the instant when the data 
returns from the periodic to the random mode, the equalizer 
tap settings will be far from their optimum (for random data) 
values, and the distortion at the equalizer output could be 
larger than the channel distortion. Hence it is desirable that 
if there are any periodic sequences in the data, it is to be 
mapped to a sequence of larger periodicity so that its spectrum 
is denser and hence equalizer is kept 'trained*. 

1.5 Timing Recovery from Zero Crossings f~14~l 

Timing recovery is one of the most critical function 
performed by a synchronous modem in data communication system. 
There are various schemes. But most modems derive their timing 
frequency from amplitude transitions in the received signal. 

One such scheme is shown in Figure 1.9. 

The clock recovery network consists of a voltage 
controlled oscillator (VCO) and a phase comparator consisting 
of the phase comparison logic and transistor controlled current 
switches. When the input s(t) =0, the one shot multivibrator 
outputs a pulse of duration 0.5 T, where T is one element 
duration. This triggers phase comparison logic circuit. The 
correction or error signal comes out of the phase comparator in 
the form of current I c . This controls the charging and dischar- 
ging of the capacitor and the voltage across the capacitor 
controls the VCO, generating the clock signal. 

From the timing diagram it is clear that at time t^, 
the clock signal is in phase and between time t^ and t£ the 
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clock signal drifts by a small amount. As s(t) goes negative 
at tj, the one shot is triggered and it produces a pulse of 
duration 0.5 T. The phase comparison logic generates 1 two equal 
width pulses QC and QC, and the current I has a waveform with 
equal width, equal amplitude positive and negative portions. 

The net change in the charge across the capacitor is zero, the 
Yco control voltage remains constant, and no adjustment is made 
on the rate and phase of the clock signal. 

As the clock signal drifts out of phase, the phase 
comparison operation at time t 2 results in a current pulse I c 
with a more positive component. Now there is a change in the 
capacitor charge and hence a change in the VCO control voltage. 
This change in the VCO control voltage results in a correction 
of the clock phase. In the figure, the clock phase is shown - 
corrected at t^, before the next phase comparison operation is 
initiated. 

One point to be noted here is that the phase comparison 
and correction operation is initiated by zero crossings in the 
received waveform s(t). Hence the performance of this timing 
wave regenerator will degrade considerably if the signal stays 
at a constant level for long periods of time. In other words, 
such an arrangement cannot handle all possible sequences of 
data. We say such a system is not "transparent" for certain 
data sequences. Thus it is desirable to map a source sequence 
on to a line sequence which has many closely spaced transitions. 
Also during idle period there should be some line sequence 
transmitted which has closely spaced transition. 
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1.6 Precession and Its Effect on the Timing Recovery 5 
Consider a phase modulated signal s(t), given by 
-SQL, 


s(t) = S g(t - nT) cos(to t + 8 ) 
— * c n 

n=- oo 


(1.16.) 


whose idle code is 8^ = 0 for all n. In the previous section 


we have obtained 

S(w) = 


27T 2 ^ 


NT ^ 


A(k ) 


k=0 


y -1 G(kQ+ nia J 6( u - m -ka-nu ) 
s c \ s 

ca >0 (l.~14) 


n 


For convenience/ we will write Equation (1.16) as 
op , 

s(t ) = ’ |a ( g(t - nT) cos to t 

n=~ cd - n J 


(1.17) 


during idle period; where a^ is a sequence of all 1 throughout, 
This can be incorporated in Equation (1.14) by letting A(k £2 ) = 
and N = 1 . Hence 


kO 

S(to) 


2?t 

T 


n 


C(nu s ) 6 ( to - to c - n q,) 


(1.18) 


For an excess band width of less than 100%/ G(nw ) will be 

S 

zero for n ^ 0. Hence 


S(») = G(0) 6( to - to c ) 


(1.19) 


It is clear that S(<o) contains a single tone and hence it is 
not possible to obtain a timing tone during idle condition, 
refer Figure 1.10(a). 

Now let us consider how this situation can be avoided 
by employing differential phase modulation with precession. In 
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a differential phase modulation, the phase of the carrier during 

A, 

the transmission of the n message symbol is given by 


© 

n 


© 


n-1 


+ 0 


n 


W M 

where 0 n is any one of the ^ equally spaced angles; — is the 

alphabet size. The advancing of the phase of the carrier by 
2 

■jj— radians independently of any change in the input data is 
known as precession. With precession © n can be written as 


e 


n 


= © 


n-1 


+ 0. 


n 


2nTc 

M 


( 1 . 20 ) 


During idle period 0—0. Assuming = 0, we can write 


© 


n 


2mt 

M 


( 1 . 21 ) 


and has period M. Hence during idle period we can write 

r j© j«t 

s(t) = Re £e n g ( t - nT) e J 


( 1 . 22 ) 


J©. 


i.e. a. 


n 


e -n = e j(2nt/M) 


(1.23) 


A(k^ ) 


nk? 


a n e 


-jnk(2n/M) 


= S e 

n=0 


■j(2n/M) n(k-l) 


= 6 - 


k-1 


(1.24) 


Substituting this in Equation (1.14), we get 


s(w) = 2 . , G(nt s + S”* 6( “ ” “c “ FT ” n V ( 1 . 25 ) 

h' ! 

ws 

Thus the effect of precession is to offset the tones by , 
producing the spectrum shown in Figure 1.10(b). Hence when 
s(t) is squared, it provides a tone at 1200 Hz. 
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S(f ) 



f = 1800 Hz 
c 



1200 


/ \ T 
Channel 

frequency 
response 


— = 1200 syrribols/sec 

100% excess band 
width 


1800 


2400 3000 


Figure 1.10(a) Spectrum Without Precession- 


— >f (Hz) 





Figure 1.11(b) 


Spectrum With Precession 
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Now consider a situation of only 50% excess band width 
as shown in Figure 1.11(a), with = 0 . The spectrum with 
precession shifts the tones by i.e. by 100 Hz in this case 
as shown in Figure 1.11(b). Still there is only one inband tone. 
Hence for the systems with a very small excess band width it is 
not possible to recover timing frequency even with precession. 
Therefore for high speed data transmission r where the excess 
band width is quite small..- employing differential phase 
modulation, spectral enrichment of A(k!J ) is necessary to 
ensure timing recovery. 

Now let us see in the context of binary data communi- 
cation what is implied by saying spectral enrichment of sequence 
^a *y. The block diagram of a differentially phase shift keying 

* V 

modulator is shown in Figure 1.12 [12]. The differential 

encoding operation performed by the modulator is explained in 


Table 1.1. The encoding operation starts with an arbitrary 

r i . 

first bit, say 1 and thereafter the encoded bit stream <jd n j> is 
generated by 


d -i h 
n-1 n 


d . b 
n-1 n 


( 1 .260 


and the sequence £ a n ] is given by 



if a n = i 

if d = 0 

n 


(1.27) 


The encoding operation during idle condition (space idle condi- 
tion is assumed) is explained in Table 1.2. It is clear that 
the sequence ^a^ is periodic with period .as two element 

l n J 
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Figure 1.12 DPSK Modulator. 


Table 1.1 Differential Encoding 


Input sequence £ b n j. 

1 

1 

1 

CD 

1 

0 

0 

0 

1 

1 

Encoded sequence £ d n 

^ 1* 

1 

1 

0 

0 

1 

0 

1 

1 

1 

Encoded sequence ^ a n ' 

\ 

1 

1 

-1 

-1 

1 

-1 

1 

1 

1 

Transmitted phase 

0 

0 

0 



0 


0 

0 

0 

* Arbitrary starting 

reference 

bit. 








Table 1.2 Differential 

Encoding 

— Space 

Idle 

Condition 


Input sequence { k> n 


0 

0 

0 

.0 

0 

0 

0 

0 


Encoded sequence ^ dj 

: 1* 

0 

1 

0 

1 

0 

1 

0 

1 


Encoded sequence ^ 


-1 

1 

-1 

1 

-1 

1 

-1 

1 


Transmitted phase 

0 


0 


0 


0 


0 



* Arbitrary starting reference bit. 
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duration. From Section 1.1 it follows that the spectrum of 
a n t will be enriched if it is mapped to a sequence of larger 
period. 


1.7 Scrambler 

From Sections 1.2 and 1.4 to 1.6 it is clear that the 
data sequence should not consist of short periodic sequences 
as well as longer durations of constant amplitude. This is 
achieved by interposing a device called "Scrambler" between the 
data source and the modulator. Its purpose is to map any 
periodic source sequence to a line sequence of larger period 
and to ensure a line sequence with closely spaced transitions 
during the time when source sequence consists of longer durations 
of constant amplitude. A descrambler is inserted between the 
demodulator and the data sink, to undo what scrambler at the 
transmitter end did, to get back the original source sequence. 

The system block diagram is shown in Figure 1.13, where £b n l is 


Scrambler Modulator 



1 Demodulator Descrambler 



Figure 1.13 


Block Diagram of a Communication System with 
Scrambler. 
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the source sequence and ^1^ is the line sequence- If there 
are no errors , the detected line sequence will be same as 

|l n | and in that case |_'b n ^ will be nothing but £b^ , the trans- 
mitted source sequence. 

A method of scrambling binary data is shown in Figure 
1.14. A binary maximal length sequence (M-sequence) is generated 
using a feed back shift register and is added (modulo 2) bit by 
bit to the data string- The data is descrambled by repeating 
the procedure using an identical synchronised shift register. 

The source sequence is recovered at the descrambler if the state 
of the random sequence generators at both ends is synchronised. 
There is no way of automatic synchronisation and hence this is 
required to be achieved by external means , may be by providing 
additional channel or circuitry. 

1.8 Self Synchronising Scramblers 

Savage [14] proposed two types of scrambler equipped 
with M-sequence generators and having self synchronising property. 
Hence they are termed as self synchronising scramblers. The two 
types differ only in their "monitoring logic" used to detect and 
eliminate any sequence of undesirable period at the output of 
the scrambler. The scrambler and descrambler structure which is 
basic to both the types are shown in Figures 1.15(a) and 1.15(b) 
respectively. 

In the structures shewn in Figures 1.15(a) and 1.15(b)/ 
the data is a sequence of digits from the modular field of p 
elements, 0, 1 , . .., p-1 , where p is prime. Addition 

and multiplication is taken modulo p. The output of the 























22 


storage elements are multiplied by the tap constants (c q , c^, 
drawn from GF(p) . The tap polynomial - also called as 
characteristic polynomial - c(x) is given by 

c (x) = x m - c ,x m “'*' - ... - c (1.28) 

m— 1 o 

is a primitive polynomial over the field GF(p) . When there is 
no input, the output sequence of the basic scrambler is all zero 
if the initial state of the storage elements is all zero and for 
the non -zero initial state, the output sequence will be a re- 
sequence of period q = p ra -l . This period we will call as 
"cycle length". From the Figures 1.15(a) and 1.15(b) we can 
write 


n 

and 


b + (c .1 . + c 0 1 -+...+C 1 ) 

n m-1 n-1 m-2 n-2 o n-m 


A A A A A 

b = 1 - (c .1 ,+C „1 _ 4- ... + c 1 ) 

n n m-1 n-1 m-2 n-2 o n-m 


(1.29) 


(1.30) 


f A m, ^ "*j 

Assuming error free transmission, i 1 \ = <1 > and hence 

it n L n J 


A 

b 


n 


1 - (c ,1 -+...+ c 1 ) + 

n m-1 n-1 o n-m 


(c ..1 .,+...+ c 1 ) 

v m-1 n-1 o n-m 


= b (1.31) 

n 

From the Figures 1.15(a) and 1.15(b) it is clear that 
the scrambler-descrambler combination have the self synchroni- 
sation property. Any error introduced in the line sequence is 
felt as long as the values stored in the descrambler disagree 
with those stored at the scrambler; i.e. for m bit duration. 

Also a single error affects as many number of bits as the number 
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of non-zero tap coefficients in the characteristic polynomial. 
Hence in practice a primitive polynomial with least number of 
feed back coefficients is selected as the characteristic poly- 
nomial. 

Savage has analysed scrambler performance for periodic 
inputs. He has shown that "the basic scrambler when excited by 
a periodic sequence of period s will respond with a periodic -line 
sequence which has either period s or a period which is the 
least common multiple { LCM ) of s and q (LCM (s, q) ) where q is 
the cycle length. The period with which the scrambler responds 
is a function of the initial values stored in the scrambler 
storage elements, that is, its initial state, and there is but 
one such state (for each phase of the input sequence) for which 
the line sequence has period s . For all other initial states 
the line sequence has the larger period" . We will call that 
initial state which maintains the period of a sequence as "the 
critical state" for that sequence. 

For the binary case Savage has found lower and upper 
bounds for the number of transitions in the line sequence when 
the scrambler is not in the critical state. Also representative 
calculations are made to determine the spectrum of the 
scrambler output and it is shown that when the source is periodic, 
the output spectrum has n times as many tones as the unscram- 
bled spectrum and each tone has 1/n the energy, where n is the 
factor by which the source period is increased. 
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1 • 9 Other Aspects of Scrambler Studies 

Apart from these properties of scrambler for periodic 
input sequences, various other authors have established that in 
general the statistics of input sequences are improved by 
scrambler in the sense that the statistics at the output tends 
towards that of white noise. This is demonstrated by them 
through autocorrelation and/or power density spectrum. 

1.9.1 Autocorrelation Function 

Nakamura and Iwadare [lO] have done analysis of a 
scrambler equipped with a multilevel M-sequence generator, which 
is a direct extension of Savage * s binary scramblers . Through 
autocorrelation analysis they established that the scrambler 
provides considerable randomness to multilevel source sequences. 

Henrikson [6] showed that when the input sequence has 
certain defined randomness , the output digits are uncorrelated 
everywhere except possibly at symbol separation that are integer 
multiples of the cycle length. At those points the output 
autocorrelation is shown to be equal to higher order moments of 
the input sequence. 

1.9.2 Power Spectrum 

Petrovic [ll] obtained a general formula for the power 
spectrum of a binary signal at the scrambler output. He assumed 
that the input binary sequence is stationary and statistically 
independent. His sample calculations of power spectrum confirms 
the facts established by Henrikson and Kasai et al. 
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1.9.3 Cascaded and Parallel Scramblers 

Gitlin et al. [5] studied the input output relationship 
of scrambler to provide extension to Savage's results. They 
observed that when the input period is a multiple of the cycle 
length q, the output period is rq independent of the initial 
state , where r = 1, 2 , . ... Also they gave a complete descrip- 
tion of the output period of a cascaded scrambler as a function 
of the number of stages. It is shown that the output period 
does not necessarily double with the addition of a stage and 
that if a particular scrambler stage does not reach critical 
state/ no succeeding stage will do so. Also, it is proved that 
the parallel scrambler configuration can, via spectral cancell- 
ation/ cause the strength of the timing tone to be attenuated. 

1.9.4 Universal Scrambler 

Lee per [8] has shown that it is possible to 'whiten 1 
to any degree all the first- and second-order statistics of any 
binary source at the cost of an arbitrarily small controllable 
error rate. He has established that any binary source can be 
'whitened 1 to an arbitrary small first- and second-order probab- 
ility density imbalance 6 if (i) the source is first passed 
through the equivalent of a symmetric memory less channel with 
an arbitrarily small but non-zero error probability e and 
(ii) the scrambler contains M stages where 

M l+log 2 [( In 26 )/ln (1 - 2 e)] 

The structure proposed by him is shown in Figure 1.16. 
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Figure 1.16(a) Universal Scrambler. 


* W** 


Figure 1.16(b) Descrambler. 



Since any binary source can be 'whitened' by this method, he 
called it as "Universal Data Scrambler" . 

1.9.5 Systematic Jitter Suppression by Scrambling 

Let us again consider the timing recovery scheme shown 
in Figure 1.5. From the figure it is clear that any obliter- 
ations in the zero crossings of the received signal causes change 
in phase of the clock signal. This clock signal is used to gate 
the equalized received signal. Hence the obliterations in the 
zero crossings in effect causes phase modulation of the received 
signal. This is known as timing jitter. This timing jitter 
introduced at each repeate gets accumulated in a repeater chain 
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and may lead to crosstalk and distortion in the reconstructed 
analog signal. The sources of timing jitter may be classified 
as systematic or non systematic according to whether or not they 
are related to pulse pattern. Systematic jitter sources are 
mainly ISI, finite pulse width, and clock threshold offsets. 

Kasai et al. [7] showed that the systematic jitter in a PCM 
repeater chain can be suppressed by incorporating scrambling. 

1.9.6 Scrambler for Encryption 

Scrambling in general is used to mean encryption for 
security purposes achieved by adding bit by bit , a M-sequence to 
a binary data sequence. Analysis of Geffe [3 ] and Yaralagadda 
[15] established the fact that any M-sequence generator structure 
can be obtained by a cryptoanalyst without much difficulty. Hence 
scrambler in general is not suitable for serious encryption. 
However by incorporating a nonlinear operation having an inverse 
operation which can be performed at descrambler, suggested by 
Arazi [l] makes the task of a cryptoanalyst more difficult. 

1 . 1 0. Scope of the Work 

In this thesis we analyse binary scramblers for periodic 
input sequences. Various cases considered are, (a)) input period 
not a multiple of cycle length, (b)' input period a multiple of 
cycle length, and <c) input period equal to i cycle length. The 
analysis is based on D-transform techniques, using the properties 
of polynomials and partial fraction expansion of rational func- 
tions of polynomials over GF ( 2 )) . 
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Also we carry out the spectral analysis of Input and 
output sequences, analogous to Fourier analysis for signals in 
continuous linear time -invariant system theory, of 2-, 3-, 4- 
and 5-stage scramblers. This analysis is based on the spectral 
theory for linear sequential circuits (LSC) suggested by Sangani 

[13]. 

1.11. Organisation of the Thesis 

In Chapter 1, we discuss certain difficulties which are 
encountered in data communication systems mainly due to periodic 
data sequences and data pattern and how this has led to the 
evolution of self synchronising scramblers. Also we discuss 
briefly the work carried out by various authors on digital data 
scramblers . 

Xn Chapter 2, we briefly enumerate certain definitions 
and properties of Galois fields and of polynomials over these 
fields. We discuss the notion of formal derivative in Galois 
fields and its application in obtaining partial fraction expan- 
sion of rational functions of polynomials over these fields. At 
the end of the chapter we state and prove an important Lemma 
which is made use of in Chapter 3 . 

In Chapter 3, we analyse the scrambler performance for 
periodic inputs through D-transform, making use of properties of 
M-sequence generater and partial fraction expansion of rational 
functions of polynomians. Through specific illustrations we 
show how a M-sequence given as input to the scrambler, loses its 
derivative properties as for as the number of transitions and 
the cyclic autocorrelation are concerned. We show that under 
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certain conditions , for a M-sequence generated by the structure 
other than that of scrambler but of sane number of stages, the 
scrambler acts just as a delay element. Also we show that the 
partial fraction expansion used throughout in this chapter can 
have two forms, one complementary to the other. At the end of 
the chapter we establish that the response of the scrambler for 
a sequence in critical state can be obtained from the responses 
due to its constituent sequences when the scrambler is in 
critical state, and that the critical state for the sequence 
will be the same linear combination of critical states for 
constituent sequences. This important result helps us to carry 
out scrambler analysis in terms of results pertaining only to a 
set of basis functions. 

In Chapter 4, we obtain basis functions for the vector 
space of periodic sequences of period equal to the cycle length 
of the scrambler and whose period remains invariant by scramb- 
ling. We give a method of obtaining Fourier series expansion 
for a periodic binary sequence over appropriate extension fields. 
We briefly discuss the existence of eigensequences and eigen- 
values, over an extension field, for LSC defined over finite 
fields . Based on this we carry out spectral analysis of periodic 

input and output sequences of 2-, 3-, 4- and 5-stage binary 

2 

scramblers with eigensequences over extension fields GF{2 ), 

3 4 5 

GF (2 )'Q GF(2 )' and GF(2 ) in each case. Detailed analysis of 
5-stage scrambler is given because the two structures in the 
case of 3- and 4-stage scramblers are reciprocal to each other 
and are quite likely to give results which probably cannot be 
generalised. 
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The last chapter contains conclusions and suggestions 
for further worlc on various topics discussed in this thesis. 



CHAPTER 2 


MATHEMATICAL PRELIMINARY 


In this chapter, first we briefly enumerate certain 

definitions and properties of Galois fields and of polynomials 

over these fields, in Sections 2.1 and 2.2. In Section 2.3 

we explain- two methods of constructing an extension field 
rn 

GF(p ) starting with a ground field GF(p) and an irreducible 
tin 

polynomial of m degree over GF(p). Partial fraction expansion 
of polynomial fractions over any field is discussed in detail 
in Section 2.4. An alternative method — similar to the one 
generally employed for polynomial fractions over the field of 
real and complex numbers — for polynomial fractions over Galois 
fields is illustrated in Section 2.5. Material for this 
chapter are mainly taken from Booth [2] and Gill [4]. At the 
end we state and prove a Lemma which is used in Chapter 3. 

2.1 Galois Field 

A field is a set of elements with two operations say, 
addition and multiplication between any two elements defined 
such that the set of elements form an abelian group under 
addition and the set of non-r^ro elements form a multiplicative 
abelian group. For example if p is a prime, the set of 
integers 0, l,...,p-l with addition and multiplication of any 
two elements taken modulo p, form a field. Such a field is 
represented by GF(p) and is called "Galois Field" of p 


elements. 
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Example 2.1: The set (0,1) with addition and multiplication 

taken modulo 2 constitute GF(2) . The addition and multipli- 
cation tables are given-.in Table 2.1(a) and Table 2.1(b). 


Table 2.1(a) Addition Table 
for GF(2) 

+ 01 

0 0 1 

110 


Table 2.1(b) Multiplication 
Table for GF(2) 

x 0 1 

0 0 0 

10 1 


Example 2.2: The set (0,1,2) with addition of multiplication 
taken modulo 3 constitute GF(3), The addition and multiplica- 
tion tables are given in Table 2.2(a) and Table 2.2(b). 


Table 2.2(a) Addition Table 
for GF ( 3 ) 

+ 012 
0 0 12 

112 0 
2 2 0 1 


Table 2.2(b) Multiplication 
Table for GF(3) 

x 0 1 2 

0 0 0 0 

10 12 

2 0 2 1 


Order , of an Element : An element p of a field is said to be 

of order k if k is the least non-zero integer such that P = 1. 
The order of any element is a divisor of the number of non -zero 
elements in the field. 

Primitive Element: If in a field, an element generates all 

the other non -zero elements by repeatedly multiplying by itself 
then that element is called "primitive element" . Its order is 
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equal to the number of non-zero elements . For example In GF(3) 
mentioned above, 2 is a primitive element. 


Characteristic of a Field: For any element in a field, if 

+ + ... + {p times) = 0, then the field is said to be of 

characteristic p. 

Extension Field GF(p m ): For every p, a prime and m an integer, 

there exist a field at (p m ) with p m elements. These are the 

m 

only finite fields existing. GF(p ) is called the extension 
field of the ground field GF(p). Any GF(p m ) is of characteristic 
p. Also any GF(p n ) CZ GF(p m ), for n, a divisor of m. 

2 

Example 2.3: The set (0,1, a, cc ) with addition and multipli- 

cation as defined in Table 2.3(a) and Table 2.3(b) constitute 
GF(2 2 ). 


Table 2.3(a) Addition Table 
for GF(2 2 ) 


+ 

0 

1 

a 

2 

a 

0 

0 

1 

a 

2 

a 

1 

1 

0 

2 

a 

a 

« 

a 

a 2 

0 

1 

2 

a 

2 

a 

a 

1 

0 


Table 2.3(b) Multiplication 

Table for GF(2 2 ) 

x 0 l a “a 2 

0 0 0 0 0 

1 0 1 a o? 

a 0 a a 2 1 

2 2 
a 0 a 1 a 


Note that GF (2 ) is contained in GF(2 ) 


2.2 Polynomials Over a Field 

A polynomial in x over a field F is an expression of 

A 

the form f(x) = t Q + f-jX + f 2 x^ + ... + f R x , where F. 
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The degree of f(x) denoted by deg f (x) is the greatest i such 
that f^ 0. The degree of polynomial f(x) = f is zero. A 
polynomial is called 'monic ' if the coefficient of the highest 
power of x is 1. The set of all polynomials over F with 
addition and multiplication of polynomials carried out as in 
normal polynomial multiplication, with the multiplication and 
addition operations between the coefficients replaced by field 
operations, forms a ring. This ring is denoted by f[x]. 

Division Algorithm: For any a(x) and b(x) in f[x], one can 

write uniquely, 

a(x) = q(x) b(x) + r(x), 0 <: deg. r(x) < deg b(x) (2.1) 

If the only divisors of a(x) are a or aa(x) - a is an element 
from F, then a(x) is said to be irreducible. 

Greatest Common Divisor (GCD): The greatest common divisor 

of two polynomials is the monic polynomial of greatest degree 
which divides both of them. Let a Q (x) and a^x) be two poly- 
nomials then with some b o (x) and b^x) it is possible to write 

GCD(a Q (x), a 1 (x)) = b Q (x) a Q (x) +b^(x) a^x) (2..2) 

Relatively Prime Polynomials: Two polynomials are said to 

be relatively prime if their GCD is 1. 

Primitive Polynomial: Every irreducible polynomial f(x) 

indivisible by x is a factor of l-x 1 for some integer i. The 
least such positive integer i is called the exponent (e) of 
f(x). The exponent of a constant is 1. If f(x) is of degree 
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n and is irreducible then its exponent is a divisor of p n -1 . 

Any irreducible polynomial of degree n and of _ exponent p n -l is 
called a primitive polynomial. For every integer n and prime p, 
there is at least one primitive polynomial . 

4 3 2 

Example 2.4: (i) Consider f (x) = x +x +x + x+ 1 over 

4 

GF ( 2 ) . f(x) is irreducible and is of exponent 5 (^ 2 -1 = 15). 
Hence £(x) is not primitive. 

4 

(ii) Consider £(x)=x + x + 1 over GF(2). This is 

15 

irreducible and of exponent 15; i.e. £(x) divides x +1, and not 
any x^+1 for i< 15. Therefore f(x) is primitive. 

Unique Factorization Theorem: Every non constant polynomial 

can be written in the form 

f(x) = aff-^x)] 1 [f 2 (x>] 2 — [f r (x)] r ; 0 (2.3) 

where is an element from F and the f ^ (x) are distinct 

irreducible# monic, non constant polynomials over F. Except 

3c . 

for ordering, f^(x) are unique „ The [f^(x)] x are called the 
prime factors of f(x). 

Reciprocal Polynomial: The reciprocal polynomial f(x) of £(x) 

n 1 ^ 

of degree n is defined as x f(rr). The exponent of f(x) is 
same as that of f(x). The prime factors of £(x) are the 
reciprocals of the prime factors of f(x). If f(x) is irre- 
ducible ?(x) is also irreducible. If f(x) is primitive f(x) 
is also primitive. 
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3 

Example 2.5: Consider a polynomial f (x) = x + x+ 1 over 

GF(2). Its reciprocal polynomial f(x) is given by 

f (X ) — X ( — T + —+1) = 1+ X 2 + X^ 

X- 3 x 

In this case, f (x) and ?(x) are both primitive. 

Certain Other Useful Properties: Following are certain other 

useful properties of polynomials over GF(p). 

(1) If f(x) is an irreducible polynomial (other than x) of 

exponent e , then the exponent of [f(x)]- 3 is p r e where r is 

such that p r "*^< j ^p r . 

3c 3c 

(2) If [f 1 (x)] 1 / [f 2 (x)] 2 ... [f r (x)] r are the prime factors 
of f (x) (indivisible by x) and if e. is the exponent of 

K 

[f A (x)] , then the exponent e of f(x) is given by 

e = LCM(e 1 / e 2 ... e r ) (2.4) 

where LCM stands for least common multiple. 

k 

(3) The prime factors of (x-x p ) are all irreducible poly- 
nomials over GF (p) whose degrees divide k. 

m_l 

(4) The polynomial (x p “ -1) splits into p-l linear factors 
in an extension field GF(p m ) i.e. 

m . _ -j 

(x P -1) = (x - a ) (x - a ) . . • (x - a P ) (2.5) 

Hence all non-zero field elements are various roots of 
unity. 

(5) (1 - x* 1 ) divides (1 - x k ) iff h divides k y i.e. (1 - x* 1 ) 
is a factor of (1 - x ) iff h is a divisor of k. 
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(6) For any polynomial f(x) over GF(p), 

Jc k 

[£(x)] P = f(x P ) . (2.6) 

2.3 Construction of GF(p m ) 

We will now see how to construct an extension field 
GF(p m ), starting with the ground field GF(p). 

Method Is Take the set S of all polynomials in GF(p) of 
degrees less than m. Define addition between any two elements 
as normal polynomial addition with the modification that the 
coefficient of variables now get added according to field 
addition (i.e. modulo p) . Define multiplication of any two 
elements as normal polynomial multiplication/ modulo an irre- 
ducible polynomial of degree m. Then the set S forms GF(p m ) 
and satisfies all field axioms. If f(x) is the irreducible 
polynomial used, then GF(p m ) is also denoted by E[GF(p); f(x)]. 

2 

Example 2.6: Consider GF(2 ) consisting of elements 0/ 1, x 

and 1+x. The irreducible polynomial used for construction is 
2 

x + x+ 1. In this addition of two elements say x and x+1 is 
given as 

(x) + (x + 1) = 2x + 1 = 1 

since 2=0 (mod 2). Also, their multiplication is given as 

(x) (x + 1) = x 2 + x = 1 (mod x 2 + x + 1) 

The addition and multiplication tables are given in Table 
2.4(a) and Table 2.4(b). 
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Table 2.4(a) 


+ 

0 

1 

0 

0 

1 

1 

1 

0 

X 

X 

1+x 

x+l 

x+l 

X 


Addition Table 
for GF ( 2 2 ) 

x x+l 

X x+1 

X+l X 

0 1 

1 0 


Table 2.4(b) 

x 0 1 

0 0 0 

10 1 
x 0 x 

x+l 0 x+l 


Multiplicat i on 
Table for GF(2 2 ) 

x x+l 

0 0 

x x+l 

x+l 1 

1 X 


Note that GF(2 2 ) of Example 
mapping 0 — 0, 1 — > 1 , a x 


2.3 is isomorphic to this with the 
2 

and a — + 1 . 


Method II 13 : Alternatively/ we can describe an algorithm 

for generating elements of GF(p m ) from GF(p) as follows. Pick an 
irreducible polynomial p(x) of degree m and coefficients from 
GF(p). Introduce a new symbol a, and assume p(a) = 0. Then 7 
0/ 1/ a, a 2 .../ a? ~ 2 will be a set of p m elements of GF(p m ) 
such that 



(2.7) 

i j (i+i) I mod p^ -1 

a a ~ a * 

Example 2.7: Consider the ground field GF(2). Let us new 

2 o 

construct e[(GF(2); x §■ x + l] = GF(2 ). Let a be an element 
in GF(2 2 ) 

.'. p( a ) = a?+a + l = 0 

i.e. a = a + 1 
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o 

Hence the set (0, 1 , a, a ) , with addition and multiplication 

2 

as given in Table 2.5(a) and Table 2.5(b), constitute GF(2 ). 


Table 2.5(a) Addition Table 


Table 2.5(b) Multiplication 


+ 

0 

1 

a 

z 

a 


0 

0 

1 

a 

Z 

a 


for GF ( 2 ) 


1 

1 

0 

** 

Z 

a 

a 


a 

a 

2 


a 

0 

1 


a 

4 

a 

a 

1 

0 


x 

0 

1 

a 

z 

a 


0 

0 

0 

0 

0 


1 

0 

1 

a 

4 

a 


Table for GF(2 ) 
2 


a 

0 

a 

2 

a 

1 


a 

0 

2 

a 

1 

a 


Addition and multiplication tables for GF(2 ), GF ( 2 * ) and GF(2 ) 
are given in Appendix A. 


2.4 Partial Fraction Expansion of Polynomial Fractions f 4 ' 

Let us now see hew partial fraction expansion of poly- 
nomial fraction over any field can be obtained based on division 
algorithm. Consider a ratio of two polynomials. 


q(x) _ q(x) 

r (x )' b^TxT t >2 (x) 

where b-^(x) and b 2 (x) are relatively prime, 
we can write with some a^(x) and a 2 (x). 


( 2 . 8 ) 

From Equation (2.2) 


a 1 (x) b^ (x ) + a 2 (x) b 2 (x) 


alsl 


Multiplying both sides by *r^y , we get 

a., (x) a n (x) 

q(x) + r 

>'2 v -°-/ 


ifey = q(x) + 571^) q(x) 


(2.9) 


( 2 . 10 ) 





rn general if r(x) = (x) b 2 (x) 


b r (x) where every b(x) 


and bj(x) (i ^ j) are relatively prime, we can 
q(x) _ q^x) q 2 (x) q r (x) 

*— • y v -f g » " — r x 11 t * * * -r 


r (x) ~ B^TxT + b 2 (x) 


FlxT 

r 


( 2 . 11 ) 


Now consider the ratio of polynomials 


all*! = a* (x). 
r ' (x) tb(x)] k 


( 2 . 12 ) 


using the division algorithm, the following equations can be 
successively constructed. 


q ' (x) = f Q (x) b(x) + r Q (x) 
f (x) = f, (x) b(x) + r, (x) 

O 1 X 

f^x) = f 2 (x) b (x) + r 2^ 


( 2 . 13 ) 


f k-2 (x) = f k-l ^ + r k-i^ x ^ 

where deg r^(x) < deg b(x) , i=0, l,...,k~l. 

By successive substitution we obtain 
q ' (x) = f 0 (x ) b(x) + r Q (x) 

= f 1 (x) [b(x)] 2 + r x (x) b(x) + r Q (x) 

- f 2 (x ) [b(x)j 3 + r 2 (x) [b(x)] 2 +r 1 (x) b(x) + r Q (x) 


= f M [b(x)f + r k _ x (x) [Mx)]*’ 1 + 
r x (x) I>(x) + r Q (x) 


( 2 . 14 ) 


• * * 


4 - 
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Dividing throughout by r'*(x) / we obtain 


qjjx) 
r v (x) 


- f k _!(x) + + 


(x) r k-2 (x) 


[b(x)]‘ 


r o (x) 

[b(x)] 1 


(2.15) 


where deg r^(x) < deg b(x), i = 0 , 1 , k-1. 


Prom Equation (2.11) and Equation (2.15) we can state 

5 

the following. If is a ratio of two polynomials where 


r(x) 
k. 




r(x) = [p 1 (x)] [p 2 (x)] ... [p r (x)] 


(2.16) 


then, can always be expanded in partial fractions in the 

form 


q(x) 

rTxT 


r i 

p(x) + Y 2 H . 


q. .(x) 




1=1 j=l [p, (x)] J 


(2.17) 


where p(x) is a polynanial and deg q^^(x) < deg p^(x). 


Example 2.8: 
given by 


Consider the ratio of two polynomials over GF(2) 


F (x) = 


5 2 

x -f x + x 

x 6 + 1 


We can write x 6 + 1 = (x 3 + l) 2 = (x + l) 2 (x 2 + x+ l) 2 

2 2 2 

where (x + 1 ) and (x + x+ 1) are relatively prime. Hence 
from Equation ( 2 . 2 ) , 


= a^(x) (x + l) 2 + a 2 (x) (x 2 + x+ l) 2 


1 
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By inspection, a^(x) = x and ^(x) = 1. This can be obtained 
by using Euclidean algorithm also, i.e. 

1 = x 2 (x + l) 2 + (x 2 + x + l) 2 


Multiplying by F(x), we can write 


F(x) = 


5 2 

X + X -4- X 
X z* + 1 


7 4 3 5 2 

X 4- X + X X + X -4- X v 

O " " f ' M " n " * '"” 1J Tr " nn "" ' " ' _ M f - \ / 

(x + X + 1) (x + 1) 


Applying division algorithm, for the first ratio 

x 7 + x^ + x 3 = (1 + x + x 4 + x 5 ) (1 + x+x 2 ) +1 

(1 + x + X^ + X 5 ) as (1 + X+X 3 )(1 + X +X 2 ) + X 

X 7 + x^ + x 3 = (1 + X + x 3 ) ( 1 + X + x 2 ) 2 + x ( 1 + x + x 2 ) +1 

2 2 

Dividing throughout by (1 + x + x ) , we get 


x 7 + x 4 + x 3 _ , . „ . __3 . 

1 — ^ ! V = 1 + X + X + 

(1 + X + X ) 


(1 + X+X 2 ) (1 + X+X 2 ) 2 


( 2 ) 


Now applying division algorithm for the second ratio, 

2 5 . , 2 3 4 W „ 

x+x + x = (1+x + x + x ) (1 + x) + 1 


2 3 4 

and (1+x +x +x) 


(1 + x + x°) (1 + x) 


. * . (x + x 2 + x 3 ) = (l+x+x 3 )(l+x) 2 +l 


Dividing by (1 + x) , we obtain 


2,5 , - 

*. . + . * ..■ 4 * « l + X + X 3 + 7 

(1 + x) (1 + x) 


( 3 ) 
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From Equations (1), (2) and (3) we can write 


F(x) = 


X + + X 5 

1 + 


X 


1+x+x 2 (1+x+x 2 ) 2 (1+x) 2 


2.5 Partial Fraction Expansion of Polynomial Fraction Over 
Galois Fields 

Booth [2] suggested a method of achieving partial 
fraction expansion of polynomial fractions over Galois fields. 
The method is similar to the one used for polynomial fractions 
over the field of real and complex numbers. 

Consider a polynomial fraction F(x) as 


F (x) 


.. A(x) 

x six) 


m m-1 

a x + a -X + 
m m-1 


+ a.x + a 
1 .. o 


x n + b ,x n ”"^ + ... + b.x + b 
n-1 i o 


(2.18) 


where m <1 n and the coefficients a. and b. are real numbers. 

JL 

The denominator B(x) can be factored into the form 

© © 0 

B(x) = (x - oc^) ^ (x - ^ ... (x - a r ) r (2.19) 


where the roots at are either real or complex numbers, 
partial fraction expansion of F(x) then is given by 


The 


F(x) = x 


1,1 

(x - 0^7 


X,2 


k- 

x l,e. 


(x - ai >' 


4 - • ♦ • 4 - 


(x - op 


+ . . . + 


JL^ 


r,e. 


(x - a r ) 


4 * • • * 4 - 


(x - a r ) 


where k, . - 


i,j (e,-j)i 


(e,-j) 


dx <e i- j) 


rr ( (x - a , ) 


e i A(x) 


B(x) ' 


(2.19) 


( 2 . 20 ) 


x « 


«L 


j = 1, 2, ... ,e^ 
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To do a similar thing for polynomial fractions over Galois 
field, it is required to define a special form of differentiation. 
In this context, let us see the meaning of the following opera- 
tion 


F(x) 


<3Q(x) 

dx 


GF(p n ) 


F(x) and G(x) are over GF(p n ) 


It means the following: 

1. G(x) is differentiated, as indicated, but the coefficients 
of G(x) ^are assumed to be arbitrary constants independent 

of x and GF(p n ) while the differentiation operation is taking 
place . 

2. The inside of the brackets is then reduced to simplest form 
without applying the addition or multiplication properties 
of GF(p n ). 

3. After step 2 is completed, the coefficients are reduced by 
using properties of GF(p n ). This operation is called the 
"formal derivative" of G(x).- 

It is extremely important that the operations are 

performed in this order. 


Example 2.9: As an illustration of the steps mentioned above 

consider the evaluation of the following expression. 


x__ d 
dx 


2i _,_2 


(■ 


lx 


(X + 1)' 


r) 


GF(2) 


x d , (x+1 ) 
21 dx 1 


x 2 (x+1 ) 


(x+1) 


GF(2) 


x___ d_ (l-x z ) 
21 ^ (x+1) 4 


GF(2) 
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F(x) 


Now consider a polynomial fraction F(x) given by 

. , , a x m + a . x 111 ” 1 + ... + a, x + a„ 

A(x) _ m m-1 1 o 

B<x) r . r-1 , , , 

x +b x +...+b,x+b 
r-1 1 o 


( 2 . 21 ) 


where m < r and the coefficients a^ and b^ are from the field 
GF(p n ). The denominator B(x) can be written as 


0 0 ! 
B (x) = [^(x)] 1 [p 2 (x)] 2 ... [p^Cx)] 


( 2 . 22 ) 


where the polynomials p^(x) are irreducible over GF(p ). 

Because the polynomials above can be of any degree/ the. general 
partial fraction expansion of F(x) will have the form 


F(x) = 


A (X> 


4 * • « « 4 - 


A i ^ 

1 /3-j^ 


Cli^Cx)] 


\,l (x 

Mk tx) 


(2.23) 


where A. . (x) is a polynomial such that deg A. . (x) < deg \i. (x) . 
l/l e. 1 

First assume that [p^(x)] 1 = [x - p/j x , where 
<= GF(p n ). Then Aj; . (x) will consist of a single constant 


term 
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A. . 

j 


(e ± -j) 

TeT-jV: - (e.-jfi f F(xHx - "] 

dx x 


e * 
x - 


GP(p n ) 


X = P H 


(2.24) 

where the term inside the brackets indicate the formal deri- 
vative of a Galois polynomial. 

Next assume that 


[^ ± (x)J 1 = [x k + c^x*” 1 + ... + c x x + cj 1 


(2.25) 


Then the polynomials A. ,(x) associated with [|i.(x)] J in Equation 

i/ j 3 - 

(2.21) are given by 

n (k-1 ) 


n 


A i,j (x) 


h. (h.) J 

2 + 1 _ 


(x " U)J (x - uP n ) j 


(h.) p 

+ *** + n(k-l) J 

(x-u p ) 


(2.26.) 


where u is a root of p,^(x) in the extension field E[GF(p n ), 


^(x)] and 


h 


i.e. 


P(x) 


r (e.-j) 

1 d 1 

B . 

( F (x ) (x— U ) 1 ) 



(e.-j) 1 (e.-j) 

x dx 1 

E[GF(p n ); ^(x)] 

x=u 




(2 

2.10: Let us consider the same F(x) as in Example 

5 2 

X + X + X 

5 , 2 

X + x ■ 

4- X 




X +1 

(x + l) 2 (x 2 + 

X + l) 2 


A l,l A l,2 

, A 2/l (x) 

A 2.2 (x) 


x + 1 (x+1) 2 

X 2 + X + 1 

(x 2 + X + l) 2 



defined over GF(2). 
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We can calculate 1 as. 


1/1 


^ F(x)(x + l) 2 


GF ( 2 ) 


X=1 


, 5 2 

x ± x 4- x ( ,v2 

dx 6 . ^ x+i; 

x +1 


GF ( 2 ) 


x=l 


d x 5 + x 2 + k 

dx / 2 , . 2 

(x + x + 1) 


GF ( 2 ) 


X=1 


8 6 2 ^ 
X 4- X + X + 1 


x 8 + X 4 + 1 


= 0 


X=1 


and A- 0 as. 


1/2 


F (x) (x+l)‘ 


[x=l 


* 5 + x 2 _±x (x+1) 2 


X 6 4-1 


X=1 


5 2 

X 4- X 4- X 

4 2 , 

X 4- X 4-1 


= 1 


X=1 


To find A 2 1 (x) and A 2 2 (x) it is necessary to go to 

2 2 

the extension field e[gf(2); x „+ x- + 1j i.e. GF(2 ). Using the 
Table 2.3(a) and Table 2.3(b) we can write/ 

(x 2 4- x 4- l ) = (x + a) (x 4- a 2 ) 

here u = a and 


dx 


5 

X 4* 

2 J 

X 4- X 

(x 

6 

X 

4- l 

5 

2 

— 

X 

4- X + 

X 

4 

2 2 

a 

x + 

Cfx 4- 


GF(2 2 ) 


x = a 


GF(2 2 ) 


x - a 



48 


8,2 6 2 4 2 2 

x -fax 4-gx tax + a 

8 4 2 

x + ax + or 


= a 


fx = a 


a 


Hence from Equation (2*26) 


A 1 1 (x) 


JOL 


x + a 


a 


(x + a ) 


(x z + x + 1 ) 


ax + a x = x 


is calculated as. 


x5 l * + * (X + a) 2 

(x° +1) 


X 


a 


5 


2 


X 


x + 

X 

T 


"TT 


x 

+ 

a x 

+ a 


x = a 


_ a 
a 


= 1 


Therefore from Equation (2.26), 2 (x) can ^ calcu l ated as 

1 1 


A. _(x) = 

1/2 1 (X .+ <X>‘ 

Now we can write 


, , 2.2 
(x + a ) 


[x 2 + x + l] 2 = 1 


F (x) = 


5 2 

X + x + X 

X 6 + 1 


X 


(x 2 + 1) x 2 + x + 1 (x 2 + X + l) 2 


This is same partial fraction expansion as obtained in Example 

2 . 8 . 

2.6 An Useful Result on Polynomials Over GF(p) 

We state and prove the following Lemma which will be 
used in Chapter 3. 


Lemma: The polynomial (1-x }, over GF(p) does not have any n 
degree primitive polynomial as its factor if k / aq where q = 
p n -l and a is some positive integer. 
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Proof : 

Case (i.) : k < q \ 

Let { 1 -x ) has an n degree primitive polynomial f (x ) 

V 

as its factor, i.e. f(x) divides 1-x for k< q which is a 
contradiction to the definition of primitive polynomial. Hence 
f(x) is not a factor. 

Case (ii): k> q, k ^ aq 

Let 

3c 1c 

(l-x^) = [f 1 (x)] 1 [f 2 (x)] 2 ... [f^(x)] -* (2.28) 

and assume that some of the f.(x)'s say f.(x), f, (x) ... f (x) 

JL JL JC T 

where i, k, . .., r j, are primitive and of degree n. Any such 

k. i 

[f^x)] has an exponent r'q where r' is such that p < k./< 
r ' 

p . Let 

k. k, k 

p(x) = [f ± (x)] 1 [f k (x)] * ... [f r (x)] r (2.29) 

Then the exponent of p(x) is given by 
e = LCM(e^/ e k , 

= lCM(rjq, r£q, . .., r^q) = r"q (2.30) 

Since p(x) is a factor of (1-x ) either k = r"q or k > r"q. If 
k > r”q, (l-x k ) should contain (l-x r ^ ) as its factor. For this 
to happen, k should be equal to sane r'“ (r"q) = aq which is a 
contradiction. Hence we conclude that (1-x ) where k ^ aq does 

i.*u 

not contain any n degree primitive polynomial as its factor. 



CHAPTER 3 


SCRAMBLER ANALYSIS - INPUT OUTPUT RELATIONS 


In Chapter 1, we have discussed that one of the main 
functions of the scrambler is to map a shorter periodic input 
sequence into an output sequence of larger periodicity. In this 
chapter we establish that the scrambler does so in most of the 
cases. From Section 3.1 to Section 3.5 we develop background 
material based on D (delay) -transform and M-sequence generator 
properties. In Section 3.6 we analyse the scrambler for its 
output periodicity when the period of the input sequence is 

(a) not a multiple of cycle length, 

(b) a multiple of cycle length and 

(c) equal to cycle length: 

of the scrambler structure. In Section 3.7 we give specific 
illustrations to bring out some interesting facts, regarding the 
nature of input sequences of period equal to cycle length and 
resulting in output sequences of same period irrespective of 
initial state of the scrambler. Also it is shown that a M- 
sequence loses its desirable properties due to scrambling, as 
for as the number of transitions and the cyclic autocorrelation 
is concerned. In Section 3.8 we discuss complementary partial 
fraction expansion of the output transform. In Section 3.9 it 
is shown that the response for linear combination of sequences 
in the critical state*, can be obtained by the same linear 


* Refer Section 1.8. 




51 


combination of the responses due to constituent sequences when 
the scrambler is in the critical state for these sequences. 
Also it is shown that the critical state for sequences can be 
obtained by the same linear combination v of the critical states 
for the constituent sequences. 


3.1 D-Transf orm 

Consider a sequence x{nT) such that x(t) = 0 for nT •<. 0 
and talcing values only at discrete time intervals . The D- 
transform o^ this sequence is defined as 


v 

■\ 

D [x(nT)] 


X(d) 



(3.1) 


where x^ = x(iT) , i = 0.1,... and x^ <~ (0,1). 

A finite sequence will have degree of d, finite. Now 
consider a sequence x(nT) to be periodic with period equal to q 
(and not any thing less than q) , i.e. 


x (nT ) — x Q ,x 1 ... Xg_^ , x o x 1 ... ... (3.2) 

Then 

D[x(nT)] = X(d) = 

= (x d° + x.d 1 + ... + x .d^Ml + d q + d 2q + ...) 
o l q— l 


_ o n q— 1 

xd +x.d + ...+x + x + . . . 

o 1 q— 1 o 




X(d) 


x d° + x,d^ + ... + x 1 d^”^ 

i q- 1 - 


1 + d 


CENTRA 1 .., : UBRAm 


X' (d) 


4- cL ’ 


(3.3) 


where deg X* (d) ^q-1. 
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The inverse D-transform of any X(d) is the sequence 
obtained by writing down the coefficients of d 1 / i = 0,1,2,..., 
in order, i.e. 

D ^[X(d)] = x(nT) = x q x^ *♦* (3.4) 

Example 3.1: Let x(nT) = 010110,010110,... Then 

X(d) = d + d 3 + d 4 + d 7 + d 9 + d 3 ^ + d 13 + d 15 + d 3 " 6 + ... 

* (d + d 3 + d 4 ) (1 + d 6 + d 12 + . . . ) 

i.e. x(d) = 3 ± d3 + 

1 + d 6 

1 . A 2 ^ ^, 3 . ^ 

Example 3.2: Let X(d) = * ■ — - • » - 

1 + d 

By long division we can write 

2 1 4 7 Q 10 11 

x(d) = 1 + d + d + d + d + d + d + d + ... 

i.e. 

X(d) = (1 + d 2 + d 3 + d 4 ) (1 + d 7 + d 14 + d 21 + .. .) 

Hence we get 

x(nT) = 1011100,1011100,... 


3.2 Re present at i chi of a Sequence by Polynomial 

A sequence f(nT) can be expressed as a polynomial in x, 
where x is an indeterminate, as follows. Let 


- f o f l 


f (nT) 
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Then# this can be expressed by a polynomial as 
£(x) = f o x° + j^x 1 + ... + £ r x r + ... 

= f.x 1 (3.5) 

135 1 

From Equation (3.1) and Equation (3.5) it is clear that the D- 
transform is nothing but a polynomial representation o£ the 
sequence . Hence they are used interchangeably in subsequent 
sections . 

Consider a sequence f (nT) obtained by bit by bit 
addition of two sequences x(nT) and y(nT) i.e. if 

x(nT) ~ x q x 1 x 2 ... 

y(nT) = y oyi y 2 ... 

then f (nT) = x(nT) + y(nT) = f o £ l f 2 

where = x ± + y ± , i = 0#1#2#... 

The transform of f(nT) is given by 

F(d) = f Q d° + f 1 d 1 + ... + f r d r + ... 

= (x q + y Q )d 0 + ... + (x r + y r )d + ... 

= x d° + x, d 1 + . . . + x d r + . . . + 

O 1 r 

■ y o d° + y^ 1 + ... + y r a r + ... 

\ 


X(d) + Y(d) 


(3.6) 
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3 .3 Periodicity o£ Two Periodic Sequences Added Bit, by Bit 
Consider two periodic sequences of period s and q, 
given by 

y(nT) = ... y s „i-Y 0 Fi ••• 

x(nT) = xx ... . xx, ... 

o 1 q-1 o 1 

Then, bit by bit addition of these two sequences can be struc 
turally represented by shift registers as in Figure 3.1. 



Let £<nT) be the resultant sequence. Obviously the periodicity 
of f(nT) depends on the periodicity of states of the circulating 
registers. Circulating Register I will be in the same state at 

any instant t = ks. Also Register II will be so for any t = k.'q. 

* 

Hence f(nT) starts repeating from the instant which satisfies 
both of these conditions. This instant t obviously is the LCM 
of s and q. Hence periodicity of f (nT) is the LCM of the 
periods of x(nT) and y(nT). 
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3.4 M-Sequence Generator f9 1 

Consider a M-sequence generator as shown in Figure 3.2. 



Figure 3.2 M-sequence Generator Structure. 

The characteristic equation of the structure is given by 

c(x) = x n + c .x 11-1 + ... + c.x 1 + c (3.7) 

n— I i o 

til 

where c^ and GF(2 ) and c(x) is a primitive polynomial of n 
degree in GF(2). 

The 2 n ~l distinct sequences of periodicity 2 n -l obtained 
n C»T. i-K'O-l 

with 2 -1 distinct non-zero er±tica -l states and a sequence of 
all zeros form a vector space of dimension n. If these sequences 
are expressed by polynomials then the least degree polynomial 
among them is called the generating polynomial g(x) and is given 
by 

g(x) = * q ± - (3.8) 

c(x) 

where q = 2 n -l and “c(x) is the reciprocal polynomial of c(x). 

The polynomial g(x) corresponds to a sequence over one period. 

If we write the sequence as an infinite sequence, then the 
polynomial corresponding to that is given by 
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g' (x) = (i + x 9 +x 2 <3 + ...) 

c(x) 

= * L ±-1 ™2 * - 1 — ( 3 , 9 ) 

c(x) x q + 1 c(x) 

3.5 Transfer Function o£ Scrambler 

Now consider the following two structures. One scrambler 
structure and another M-sequence generator with their initial 
states indicated in Figure 3.3. 



It is clear that the output sequences from the two structures 
are sane; i.e. we can say that the impulse response of the 
scrambler is nothing but a M-sequence. From the M— sequence 
structure it can be seen that the M-sequence generated corres- 
ponds to the least degree polynomial (over one period), i.e. 
g(x). Hence from Equation (3.9) we can write H(d), the transfer 
function of the scrambler as 

H(d) = -i- = (1 i-q- = <3-10) 

c(d) ’S(d) (1 + d q ) 1 + <T- 

where deg q(d) = q - n. 
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Since all M-sequences plus one all zero sequence form a n- 

rv 1 

dimensional vector space, H(d) , dH(d) , d H(d) corresponds 

to n linearly independent sequences of that vector space. 

Hence we can write the transform of any M-sequence as 

p(d> = s(d) H{d) = s(d) ■ (3.11) 

1 + d^ 

where s(d) is some polynomial of degree <C n-1 . Since (l+d g ) 
has all the irreducible polynomials except d of degrees which 
are divisors of n as factors, it is clear that q(d) will have 
all such factors except c(d) . 

3,6 Output Period of Scrambler 

Consider the scrambler structure shown in Figure 3.4. 




f y} 


Figure 3.4 Scrambler Structure.. 

Let X(d) be the input transform and H(d) be the transfer 
function. It can be shown that the system is linear. Hence 
we can write the output transform Y(d) as 


Y(d) 


X(d) H(d) 


(3.12) 





&a 


New let us study the response of the system to periodic inputs. 

Case I: Let the input be of period s where s ^ rq, r = 1,2,... 

The input transform be given by 


X(d) 


1 + d" 


deg p{d) < s-1 


(3.13) 


Then Y(d) = x(d) H(d> 


(1 + d ) (1 + d?) 


(1 + d s ) ~(d) 


(3.14) 


Let (1 + d ) = [f. (d)] 1 [p„ (d)] z ... [f (d)] 


(3.15) 


From the Lemma proved in Section 2.6, we know that ?5(d), being 

til s 

an n degree primitive polynomial, is not a factor of (1 + d ). 

Hence, we can write Y(d) by the partial fraction expansion method 


A 1 1 (d) 

Y(a > = -V5T + 


A i,*< d) 

+ + ... + 

[*!«>] 1 


A r,l (d> 

V d) 


+ . . . + 


A r k (d) 
r,k r 


0M d) l 


_ + A (d) 

c r 'c(d) 


(3.16) 


By combining all terms other than the last one, we can write 


Y(d) = 


p* (d) + A(d) 

1 + d s 'c(d) 


1 + d^ 


A (d) g(d) 


1 + d- 


deg A(d) n-1 


(3.17) 


This is the output transform,, only if the initial state is 
zero. For non-zero initial states, transform corresponding to 
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zero input response is to be added to Equation (3.17). 
from Equation (3.11) we can write Y(d) as 

Y(d) = p* i£l + A(d) q(d) + s(d) q(d) 

1 + d s 1 + d q 1 + d q 


Hence 


(3.18) 


where deg A(d) , deg s(d) 4 n-1. 

i.e. Y(d) - . q(a) CA(d) ± 5(d)] <3 . 19) 

1 + d s 1 + d q 

It can be seen that there will be one particular initial state 
for which A(d) = s(d). In that case, 

Y(d) « (3.20) 

1 + d S 

i.e. the scrambler is in a critical state for that particular 
sequence and the output period is s itself. For other initial 
states 


Y(d) » , deg q'(d)j* q-1. (3.21) 

1 + d s 1 + d q 


i.e. the output sequence is of periodicity LCM(s,q). 

From. Equation (3.20) and Equation (3.21) we can state 
that if the input sequence to a scrambler is of period s which 
is not a multiple of ( 2 n -l), then the output sequence is of 
period LCM(s, 2 n -l), but for a unique initial state for which 
the output period is s itself. 


Case II: Let s = rq then. 


Y(d) = 


p(d) q(d) 

(1 + d S ) (1 + d q ) 


deg p(d) ^ rq- 1 , 
deg q(d) = q-a 


P.l±L3l£l 

(1 + d rq )(l + d q ) 


(3.22) 
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(a) If p(d)q(d) has either l4d rq or l-fd q as a factor then y(d) 
is given by either 


Y(d) 

_ P 1 (d) 

(3.23) 


1 + <3 q 


Y(d) 

- P" (d) 

(3.24) 


1 + d rq 


i.e. the output period is either q or rq. 

(b) Let us write Y(d) as. 


Y(d) * — — (3.25) 

(1 + d rq ) c(d) 

Since (1 + d q ) is a factor of (1 + a rq ) , *S(d) is definitely 
a factor of l4d rq . If p(d) contains this factor, then 
there will not be a repeated factor c(d) in the denominator. 

(3.26) 

LCM(q, rq) = rq. 


In that case, we can write 


Y(d) = 


p'<a> + g'(d) 


1 + d' 


rq 


1 + d - 


Therefore the output period = 


i.e, 


Y(d) 


p" (d) 


(3.27) 


(c) Let us say, after cancellation of common factors if any. 


Y(d) 


p' (d) 


}c 1c }c 

fF. (d) ] 1 [F ~ (d) ] 2 ... [F (d)] r 


(3.28) 


If we denote the exponent of F. (d) by e, , then the exponent 

m i 

of [F. (d)] is given by p e^ where nu is such that 
m.-l 1 m. 

p :L <^ ^ p 1 and the exponent of the entire denominator 
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m l m o m_ 

is given by LCM(p e^, p ^e 2 , *** P e r ) * I£ there is a 
primitive polynomial (at least one) of degree n as a factor, 
then , 


m l m 2 rn _ 

LCM(p e 1 , p e 2 , . .., p e r > 


Hence Y(d) can be written as 


r 'q, r ‘ £= (1,2, . .. ) , 

(3.29) 


Y(d) = 


p" (d ) 

(1 + d r ‘ q ) 


(3.30) 


i.e. the out-put sequence is of period r*q. 

The cases discussed in (a), (b) and (c) above are all 
for zero initial state. For non-zero initial states, Y(d) can 
be written as. 


p" (d) + s (d) q(d) 

1 + d r * q 1 + d q 


(3.31) 


Hence the output period will be equal to LCM(r'q, q) , i.e. r'q 
itself . 

Hin Equation (3.28), if none of the F.(d), is a primitive 

m^ m 2 m 

polynomial of degree n, then the LCM of (p e-^, pe 2 , ..., p r e r ) 
need not be of the kind r'q. Let us prove that this cannot 
happen and hence the output period is always some r'q. 

If the exponent of the denominator in Equation (3.28) 
is some q', not a multiple of q, then 

Y(d) = - — — - T where deg a(d) ^ q' - n - 1 

1 + d q 


and from Equation (3.12) 
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X(d) 


= a(d) 1 ' _ a(d) 'g(d) 


1 + d q ’ 


1 + d 


q' 


i.e. X(d) 


a«(d) 
1 + d q 


deg a* (d) < q - 1 


(3.32) 


This contradicts our assumption that the input period s is equal 
to rq. 

We can state that when the input is of period rq, the 
output period will be some r*q, irrespective of the initial state. 

Case III: Let s = q 

Though this is covered in Case II, let -us consider 
separately the case of input period equal to q, i.e. let 

if 

f 

* 

X(d) * , deg p(d) ^ q - 1 

1 + d q 

Then Y (d) = X(d) H(d) = -1— ^ HjL 3!lAL- (3.33) 

1 + d q 'c(d) \-t-J* i+JV 

If p(d) q(d) has a factor (1 + d q ) , then the periodicity of the 
output will be q. If p(d)q(d) does not have (1 + d q ) as a 
factor, then the output period will be equal to 2q. Let us 
say p(d)q(d) has (1 + d q ) as a factor i.e. 

p(d) q (d) = r(d) (1 + d q ) 

p(d) = U ^qfdj = Z(d) r(d) (3.34) 

where, deg p(d) ^ q - 1? deg c(d) = ru Therefore, 
deg r(d) < q - (n + 1) = 2 n - n - 2 
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/ -4 \ 

Hence there are (2 l “ n ~ x ; _ 2 ) sequences of period q, which 
gives rise to output period of q itself (excluding all zeros 1 
and all ones 1 sequences ) . 

Substituting for p(d) in Equation (3.33) 


Y(d) 


c(d) r(d) 1 
1 + d q ’ c*(d) 


r(d) 


1 + d ■ 


(3.35) 


If the initial state is non -zero, then 


Y(d) 


r(d) + s(d) q(d) 
1 + d q 1 + d q 


r 1 (d) 

1 + d c 


deg r 1 (d) ^ q - 1 


(3.36) 


i.e. the output period is q, irrespective of the initial state 
for these (2^ 2 “ n “ 1 ^ - 2) sequences except when p(d> corresponds 
to a sequence with internal periodicity. For these cases if 
we write p(d) as 


p(d) = (Pq_i dq_1 + Pq-2 dq 2 + *♦• + P^ + P Q ^ (3.37) 

then it can have periodicities of period i, for all i‘s which 
are non-trivial divisors of q. The total number N of such 
sequences is given by 

(2 d - 2) (3.38) 


These sequences will give rise to period i itself when scrambler 
is in critical state, with this, we can state the following. 

For a scrambler of cycle length q where q = 2 — 1 , 
there are exactly N* number of input sequences of period q. 
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which results in output sequences of period q itself irrespective 
of initial state, where N' is given by 

N« = 2* 2 ~ n-1 } - 2 - 2_, (2 1 “ 2 > (3.39) 

ik 

Example 3.3: Consider a scrambler structure of characteristic 

3 

polynomial c(x) = x + x + 1 



Scrambler Structure for Example 3.3. 


Impulse response is, 1011100,1011100 


.*. H(d) = 


i -.2 ,3 -4 

1 + d + d 4- d 

1 + d^ 


The zero input response and their transform for various non-zero 
initial states are given in Table 3.1 

Let the input transform be given by 


X(d) 


d 2 + d 

it? 


then Y(d) = X(d) H(d) = 


d," 


d + d 1 + d + d ± d 


1 + d~ 


1 + d 


(d 2 + d + l)(d 3 + d 2 + 1) 


i.e. Y(d) 
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Table 3.1 Zero Input Response and Their Transform 


Initial 
X 1 X 2 

t 

State J 

X 3 ; 

Response 

t 

1 

; Transform 

t 

s 

I 

1_ - — - 

0 

0 

1 

1011100,10 ... 

H(d) 

0 

1 

0 

1110010,11 ... 

(1 + d)H(d) 

0 

1 

1 

0101110,01 ... 

d H(d) 

1 

0 

0 

0111001,01 ... 

(d + d 2 )H(d) 

1 

0 

1 

1100101,11 ... 

(1 + d + d 2 )H(d) 

1 

1 

0 

1001011,10 ... 

<1 + d 2 )H(d) 

1 

1 

1 

0010111,00 

d 2 H(d) 


From Section 2.5, we can write this as 


A. - (d) A (d) 

Y(d) = + L l ±— 


2 3 2 

d + d + 1 d + d + 1 


( 1 ) 


where 


A l,l (d) (d + a) 


h l h l 

4 - 


d + a 


(a 4- d + 1 ) , 


( 2 ) 


is a root of d + d + 1 = 0 in E[ GF (2); d + d + 1] , 


h. 


d 


(d 2 + d + 1) (d 3 + d 2 + 1) 


. (d + a) 


GF(2 2 ) 


d = a, 
(3) 


and A 2 ^ 


h_ h; 

1 + 5T- + 


h; 


(d + P } (d + B 2 ) (d + P 4 ) 


(d 3 + d 2 + 1) 


is a root of d 3 + d 2 +1=0 in e[gf( 2 ) ; d + d + ll 


( 4 ) 
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h « — 2 2 . (d + P ) 3 

j_(d + d + l)(d + d + 1) jGF(2 J )|d = p 


Now let us calculate .^(d) and A 2 1 (d). From Equation (3) 
and referring to Table 2.5(a) and Table 2.5(b), 


(d + a) (a + d 2 +1) GF(2 2 ) d =a 


d 4 + ad 3 + a 2 d 2 + a + a 2 d = a 


= a 


From Equation (2), 


0£- fat 

A. 1 (d) = -r— - + o 

1,1 d + a d + a 2 


(d* + d + 1) = d + 1 


3 3 2 

From Equation (5) and referring to tables of e[gf( 2 ) ; d + d + Ij 
in Appendix A 


(d 2 + d + 1) (d + p 2 ) (d + p 4 ) |GF(2 3 ) d = p 


Since d 3 + d 2 +1 = (d + p ) (d + p 2 ) (d + p 4 ) 


12 d 4 + pd 3 + P 2 d 2 + p 3 d + p 6 d = p 


Therefore from Equation (4), 

A ) — § — + — § — s- + — _ = d 2 + d + 1 

A 2,l (d > <J + e + d + p 2 d+p 4 


Hence we can write Y(d) as 
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Y(d) = 


d + 1 . d +d + 1 

T 


d 2 + d + 1 (d 3 + d 2 + 1) 


1 + d 2 + (d 2 -f d + 1)(1 + d)(l + d + d 3 ) 
1 + d .1 + d 7 


i.e. Y(d) 


o + (d + d + 1) H(d) , 

1 + d 


for zero initial state. From Table 3.1 we can see that for 
initial state 1 0 1, 

2 

Y(d) = + (d 2 + d + 1) H(d) + (d 2 + d + 1) H(d) 

1 + d^ 


i.e. the output period is 3 itself and the output sequence 
y (nT) is , 101,101,... 

For any other initial state. 


Y(d) 


1 + d 2 , s (d ) ( 1 + d 2 + d 3 + d 4 ) 

3 " i " 7 

1 + d 1 + d 


deg s(d)<C 2 


Hence output period is LCM{3,7) =21 


Example 3.4: For the same structure as in Example 3.3, let 


X(d) 


1 + d + d 3 + d 4 + d 5 + d 6 + d 9 


21 


1 + d 


(1 + d 2 + d 3 ) 2 (d 3 t d t 1) 
21 


1 + d 

We have Y(d) * X(d) H(d) 


(1 + d 2 + d 3 ) 2 (d 3 -f d + 1) 


i.e. Y(d) 


1 + d 


1 + d 

df 

21 


21 


3 2 

d+d+1 


(1 + d 2 + d 3 ) (d 3 + d + 1) 
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_ _ 1 

(1 + d)(d 2 + d + l)(d 6 + d 4 + d 2 + d + l)(d 6 + d 5 + d 4 + d 2 + 1) 

Here/ (X + d) is of exponent 1 
2 

d + d + 1 is of exponent 3 
6 4 2 

d + d + d + d + 1 is of exponent 21 
d + d + d + a + 1 is of exponent 21 

Hence exponent of denominator is 104(1/ 3/ 21/21) = 21 and we 

can write 


Y (d) = 


(d 3 + d 2 l)(d 3 + d + 1) 


1 + d 


21 


d 6 + d 5 + d 4 + d 3 + d 2 + d + 1 


i + a 


21 


i.e. the output period is 21 and the output sequence y(nT) is, 
111111100000000000000,111 ... 


3.7 Specific Illustrations 

Based on the general analysis we will now consider the 
behaviour of scrambler with specific structure and/or specific 
input . 

(1) Consider a three stage scrambler. Prom Equation (3.38), 

N’ = 2 (8-3-1) _ 2 _ 14f since ^ t 2 X -2=0 

i 

i j 7 
i/7 

i.e. there are 14 distinct sequences of period 7 which results 
in output sequences of period 7 itself . 

Now consider a M-sequence generated with c*(x) - the 
reciprocal polynomial of scrambler structure as the characteristic 
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^ 3 2 

polynomial. Say, c(x) = x + x + 1, Then, any M-sequence 
generated by it will have the transform 


X(d) 


r‘(d)(l + d)(a 3 + d 2 + 1) 
1 + d 7 


deg r ' (d) < 2 


(3.40) 


Let this be given as input to the scrambler structure with 
characteristic polynomial c(x) = x 3 + x + 1. Then, 


Y(d) * X(d) H(d) 


r* (d)(l + d) (d^ + d 2 ± 1) 
” 1 + d 7 ~ 


(1 + d) (d + d ± 1) 
(1 + d* 7 ) 


Y(d) = 


_ r* (d) (1 ± d? 


(3.41) 


1 + d' 


For non -zero initial states 


Y(d ) = r* (d) (1 + d) + s 1 (d) (1 + d)(l + d + d~*) # 


1 + d“ 


1 + d ■ 


deg s ' (d) 4 2 


i.e. Y(d) = 


r" (d) 


1 + d 


7 ' 


deg r" (d) .< 6 


(3.42) 

(3.43) 


Therefore the output period is 7 , independent of the initial 

~k 

state. Cbviously the complements of these seven sequences also 
give rise to output period equal to 7; i.e. there are 14 such 
sequences. Hence for three stage scramblers, we can state that 
only M-sequence s generated by the reciprocal structure (with 
respect to scrambler structure) and their complements give rise 
to output sequence of period 7. 

* A sequence complement to x(nT) is given by 
x(nT) = x Q x x x 2 . . . 
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(2) In the case of 4 stage scrambler from Equation (3.32), 

N‘ = 2 (2 “ 4 " 1) - 2 - (2 1 - 2) = 2010 

i 

1 1 15 
1^15 


Now consider a M-sequence generated with c (x) - the reciprocal 
polynomial of scrambler structure - as the characteristic poly- 
nomial; say c(x) = x 4 + x^ +1. Then, from Section 3.5 any 
M-sequence generated by it will have the transform. 


X(d) = 


r * (d) (1 + d)(l + d + d 2 ) (1 + d 4 - d 2 ± d 3 + d 4 ) (1 ± d 3 +d 4 ) 


1 + d 


15 


deg r* (d) ^ 3 (3.44) 

Let this be given as input to the scrambler structure with 
characteristic polynomial c(x) = x 4 + x + 1 = .0. Prom Section 
3.5 the transfer function H(d) is given by 


H(d) 


_ (1 + d)(l ± d ± d 2 ) (1 ± d+ d 2 ± d 3 ± d 4 )(l jd ± d 4 ) 


1 + d 


15 


Y(d) 


r' (d) (1 + d) (1 + d + d 2 ) (1 ± d + d 2 ± d 3 - 4 - d 4 ) 


1 + d 


15 


For non-zero initial states. 


Y(d) = 


r ' (d) (1 + d) (1 + d + d 2 )(l + d + d 2 + d 3 + d 4 ) 


1 + d‘ 


15 


s' (d)(l ± d)(l + d + d 2 )(l ±d± d 2 -f d 3 + d 4 )(l ± d + d 4 ) 


1 + d' 


IT 


deg s 1 (d) <T 3 (3.45) 


- , deg r" (d) ^ 14 


1 + d 


(3.46,) 
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i.e. the output period is 15 independent of the initial state. 
There are 15 such M-sequences. Obviously, sequences complement 
to these sequences also give rise to output sequences of period 
15, irrespective of the initial state. 

(3) New consider a scrambler structure of n stages in general 
with a characteristic polynomial c(x). Then, the transfer 
function H(d) is given by 

H(d) = — i— (3.47) 

S(d) 

If the output of any other M-sequence structure of degree n with 
characteristic polynomial c'(x) is given as input to the scrambler, 
then from Section 3.5 we can write the input transform as, 

X(d) = , deg s(d)^ n - 1 <3.48) 

ZXd) 

Therefore output transform is given by 

Y(d) = (3.49) 

^'(d) £(d) 

Here 'c' (d) and 'c(d) are two different primitive polynomials of 
degree n. Hence Y(d) can be written as 

y (d) = , q - 2n ^ deg s’ (d) 4 q - n - 1 (3.50) 

1 + d q 

For non-zero initial states, from Section 3.5 we can write 

Y(d) - s . 'J d- L. + t q - n <C deg r (d) ^ q - 1 (3.51) 

1 + d q 1 + a 1 

i.e. Y(d) = S " ■ — • 7 - , q - 2n ^ deg - 1 (3.52) 

1 + d q 



The output period is q independent of initial state. It follows 
that this is so when the input is any sequence complement to 
these input sequences. Therefore we can state that for a 
scrambler of M stages with characteristic polynomial c(x). If 
the sequences generated by any M-sequence generator of the same 
number of stages and characteristic polynomial other than c(x) 
or the complement of these sequences are given as input, then 
the output period will be equal to q, the cycle length. 

The observation made with respect to 3 and 4 stage 
scramblers respectively in Illustration 2 and Illustration 4, 
is covered in this general statement. 

(4) We have seen that there is no advantage in scrambling a 
M -sequence of period q, as for as its periodicity is 
concerned, except when the sequence is one of the sequences 
generated by the same structure as that of the scrambler. 

The next question is whether other desirable charact eristics 
of a M-sequence, say number of transitions and cyclic 
autocorrelation function (CACF) are still retained at the 
output. 

Consider the three stage scrambler with characteristic 
3 

polynomial c(x) = x + x + 1, with input, a M— sequence generated 
by the structure with characteristic polynomial £?(x) = x + x + 1. 
Then Equation (3.42) gives the output transform Y(d) for non- 
zero initial states as 

r ' (d) (1 + d) , s' (d) (1 + d)(l + d + d 3 ) 

7 + 7 

1 + d 1 + d 


Y(d) 
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For a given r'(d) there can be solutions to s * (d) such that the 
numerator say, p' (d) is of the kind 

p*(a) * d 1 + d 1 ® 1 + ... + d 1 ® 10 (3.54) 

where 1 G (0,1 ,...6), im G (0,1,... 5) and © is addition, modulo 7. 
In that case,Y(d) corresponds to an output sequence which has 
only two transitions over one period. All such combinations of 
r'(d) and s'(d) and p 1 (d) resulting in these cases are listed in 
Table 3.2. It can be seen that there are three classes of output 
sequences having minimum number of transitions, i.e. 2. Obvi- 
ously the CACF of all sequences In a class is same. The CACFs 
for these three classes are shown in Figure 3.5. It is seen 
that the ideal CACF of a M-sequence is lost due to scrambling. 

It is verified exhaustively in the case of 3, 4 and 5 
stage scramblers that the ideal CACF of a M-sequence in almost 
all the cases, is lost due to scrambling. In most of the cases 
the number of transitions in the output sequence is less than 
that in the input M-sequence. Only in a very few cases they are 
more than that in the input M-sequence. The minimum number of 
transitions observed for four stage scramblers is 4. 

(5) Now let us consider another interesting question, i.e. 

whether it is possible to get a M-sequence at the output 
of a scrambler in a situation discussed in the previous 
illustration. Let c(x) be the characteristic polynomial 
of a n stage scrambler. Consider a M-sequence generated 
by a n stage M-sequence generator with the characteristic 
polynomial c ' (x) , as the input. Proceeding exactly the 
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Table 3.2 Cases of Minimum Number of Transitions 


r' (d) 

i 

i 

I S 1 

* 

~J 

1 

(d) i 

1 

P'(d) 

1 

1 

d 2 


2 3 4 

d + d + d + d 

1 + d + d 2 + d 4 + d 5 + d 6 


d 


d 2 4 d 3 + d 4 + d 5 

a 

1 4 

d 2 

1 + d + d 2 + d 3 + d 5 + d & 

2 

a 2 ' 


d + d + d 4 d 

d 

1 + 

d + d 2 

1 + d + d 2 + d 3 + d 4 + d 6 


1 


d 3 4d 4 

1 + d 

d 


2 3 4 5 

1 4 d 4 d 4 d 4 d 4 d 


CM 


1 4 d 4 4 d 3 + d^ 

o 

a 2 


1 4 d 4 d 3 4 d 4 4 d 5 4 d 6 

1 + d 

d 2 - 

f d 

1 4 d^ 


d 


4 5 

d 4 d 

d + d 2 

CM 


d 4 d 2 4 d 3 4 d 4 4 d 5 4 d 6 


l + 

a 2 

1 4 d 4 d 2 4 d 3 4 d 5 4 d 6 


1 ™ 


1 4 d 2 4 d 3 4 d 4 4 d 5 4 d 6 

2 

1 4* d 4* d L. 

l + 

d 2 

^ 6 
d 4* d 


d + 

d 2 

1 4 d 4 d 2 4 d 6 
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same way as we did for obtaining Equation (3.44), the 
output transform for any initial state can be written as 


Y(d) 


r(d) p(d) + s(d) p(d) c(d) 


1 + d ■ 


4 


= Hill 1 r <3) + s(d) c(d)l 
1 + d q 


where p(d) 


1 +d q 
c(d) c'(d) 


deg r(d), deg s(d)^ n-1 

(3.55) 

(3.56) 


Y(d) can c orre s pond to a M— sequence of length q only under 
the following conditions. 


(1) r(d) + s(d) c(d) = r'(d) c(d), deg r ‘ (d) <: n-1 

(2) r(d) + s(d) c(d) = r"(d) c‘(d), deg r"(d)< n-1 


It is obvious that the first condition can be satisfied 
only when r(d) = 0 and in that case r ’ (d) = s(d) and the output 
is nothing but the zero input response of a scrambler, a M- 
sequence. This is a trivial solution. There can be solutions 
to s(d) for a given r(d) such that the condition 2 is satisfied 
and in that case it is evident that the scrambler acts just as 
a delay element. Such cases for 3, 4 and 5 stage scramblers 
are listed in Table 3.3. In all cases there is only one 
solution to s (d) . 


3.8 Complementary Partial Fraction Expansion 

Let us now see whether the partial fraction expansion 
of the output transform we have used in this chapter is unique. 
After cancellation of terms let Y(d) be given by 
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Y(d) = 


A(d) 

]<r ^ 

[p x <a>] 1 [f (d)] 2 . 


• P? «>] 


(3-57) 


The partial fraction expansion of Y(d) be 


Y(d) 


A Ll (d) 


A r,l< d) 

P r (d) 


A l,k/ d) 

+ i — X" + ••• + 

[F-^d)] x 

A r,k (d) 

* + k 

[M d >3 r 


(3.58) 


where deg A. .(d) < deg F . (d) 

1 / J X 


jv . 

If the exponent of (l\(d)] 1 is denoted by e^ ^ then we can 


write Y(d) as 


Y(d) = 


B i.i <d) 


1 + d 


4 * ♦ * 4 - 


B X 

1 + d el ' k l 


4- . * • 


B r.l (d) 


1 + d 


+ ... 


B r,k r (d > 

1 + d Sr ' k r 


deg B 1/j W) e i, j - 1 


(3.59) 


Let R(d) be another polynomial fraction given by 


R(d) 


*1,1® 


l + d 


i I d ei -*i 


4 * . # . 4 


B r.l (d) 


1 + d 


4 - * * * 4 


B r,k r (d) 

1 + d>' k r 


(3.60) 


where deg- j(d) ^ e j[ / j " ^ an d A(d) is defined by 


A(d) 


.& + a.d ■+■ ... + a i d 

o i- CT ” J - 


if A(d) = a Q + a x d + ... + a -d- 


(3.61) 



79 


Their sum is given by 


A(d) + A(d) = (a o + a Q ) + (a 1 + a 1 )d + — + (a, + a 1 )d q ” 1 


■ 1 + d + ... + d q “ 1 

Adding Equation (3.59) and Equation (3.60) 


(3.62) 


Y(d) + R(d) 




4 . » * 4 - 


fci^fd) + B 1/k ^ d) ] 


1 + d 


'1,1 


'1/k. 


1 + d 


[B (d) +5 .(d)] 

+ ... + — — 

e r 1 
1 + d r/i 


4 • • * 4 


[Br^d) + B^W)] 


'r,k 


1 + d 


= TTii + tTa + •••' 3 “ 53 k j ti,nes 


(3.63) 


3=1 


Now if j is even, then 


Y(d) + R(d) = 0, i.e. Y(d) = R(d). Therefore 


Y(d) 




+ + 


B 1 v (d) 


B x (d) 

+ 4 - - — £-£-4 

« « « • » 


1 + d 


' 1,1 


*l,k. 


1 + d. 


1 + d 


'r,l 


+ . . . + 


5 r,* r (d) 


(3.64) 


'r,k 


1 + d 


is also a solution if j is even. Now consider Equation (3.17) 
i.e. 


p 1 + a.m 


y(d) „ Sim- + iWLaMl 

1 + d 1 + d y 

where deg A(d) ^ n-1, deg q(d) ^ q-1, 


1 + d s 1 + d q 
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From Equation (3.64) it follows that Y(d) can also be written as 


Y(d) = + jlM 


1+d 1 + d 


— , deg q' (d) ^ q - 1 


( 3 . 65 ) 


For non-zero initial state ,Y(d) is given by 


Y(d) 


+ BiZUim , deg s(d) £ n-1 (3. 


i + d i + d 


66 ) 


1 + d ■ 


In this case , for a particular initial state 


Q 1 (d) + s (d) q(d) _ 1 

1 + d q 1 + d q 1+d 

i.e. it corresponds to a sequence of all ones. Hence 


(3.67) 


Y(d) 


P' (d) 

1 + d S 


+ 


1 

1+d 


£ L^l 4. 2 + d + d^ + 

1 + d s 


4* . • 


i. e . Y(d) - 

1 + d s 


(3.68) 


same as obtained in Equation (3.20). 

3.9 Critical State and Output of a Composite Sequence 

Consider a sequence which is a linear combination of 
sequences. The critical state for each constituent sequence and 
the response of the scrambler in that state are known. Let us 
now investigate whether the response for the composite sequence 
can be obtained from the responses for constitutent sequences. 

Let x 1 (nT) / x 2 (nT) / ..., x r (nT) be sequences of period 
s (s ^ rq) given as input to a scrambler of cycle length q. 

From Equation (3.19), for each input the output transform for 
non-zero initial states can be written as 
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V d) 


Pj(d) q(d) A, (d) + s. (d) 
= — + — 3- 


1 +d k 


1 + d ■ 


If tills composite signal given fay 


(3.69) 


x(nT) = x 1 (nT) + x 2 (nT) + ... + x r (nT) 

is given as input/ then the output transform Y(d) for zero 
initial state is given fay 


Y(d) * Y£(d) + Y£(d) + ... + Y^(d) 


— g“ [Pi (d) + p'(d) + ... + p‘ (d) ] 

(l+d 3 ) 1 2 r 


+ “- H(d) a [A. (d) + A 0 (d) + ... + A ( d ) ] 

(1 + d q ) 1 2 r 


— Pi<d> ^ a ± 

(1 + d S ) i=l 1 (1 + d q ) i=l 1 


(d) 


For non-zero initial states we have 


Y(d) = 


(1 + d ) 1=1 


t p; + if: a (d> + s(d>] 


(1 + d' 1 ) i=l 


(3.71) 


From Equation (3.71) it follows that s(d) will correspond to 
the critical state for the composite sequence x(nT) if 

s(d) = j s. (d) (3.72) 

i=l 1 

where s^(d) corresponds to the critical state for x^(nT). The 
output transform for each x ± (nT) alone when the scrambler is in 
critical state is given by 
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P! (d) 

Y ± (d) = -1 (3.73) 

1 + d s 


Also the output transform for the composite sequence x(nT), 
when the scrambler is in critical state , is given by 

Y(d) = -- 2Z Pj(a> (3.74) 

(i + d s ) fm. 1 

From Equations (3.72), (3.73) and (3.74), we can state the 
following. 

If the input sequence x(nT) of period s (^ rq) can be 
expressed as, 

x 


(nT) 


r 

= z: 

fir 


a.x. (nT) 
r i 


a ± <c (0,1) 


(3. 75) 


then, in the critical state X c , the output sequence is given by 


Y(nT) 


2l a i^i^nT) 

i=l 


(3.75) 


where y, (nT) is the response for x. (nT) when the scrambler is 

JL 

in critical state X . and the critical state X is given by 

C 4f 1 


X. 



a. x . 
a. c,i 


(3.77) 


Example 3.5: Consider the three stage scrambler with charac- 

teristic polynomial c(x) = x^ + x + 1 . The critical states and 
outputs for four input sequences of period 5 are listed in the 


Table 3.4. 



Table 3.4 


Critical States and Responses 


Serial 
No. (i) 

, 

Input Sequences 
x ± (nT) 

Critical 

X o i 

X X. 

o 1 

States 

X 2 

Output Sequences 
(nT) 

1 

10000,10... 

1 

0 

0 

11001,11. . . 

2 

01000,01... 

0 

0 

1 

11100,11.. . 

3 

00010,00... 

1 

1 

1 

00111,00... 

4 

00001,00.. . 

1 

1 

0 

10011,10... 

5 

11011,11... 

1 

0 

0 

10001,10. . . 
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Note that x 5 (nT) is given by 


x,- (nT) = ]>*~ x, (nT) 

x~i ^ 


y 5 (nT) - ^ 


> , y 4 (nT) and 

T=f 1 


x 


4 , 

: - 5 = k **•*■ 


vrhich is true as verified from the Table 3.4, 



CHAPTER 4 


SPECTRAL THEORY OF LSC AND ITS APPLICATION TO SCRAMBLERS 

Un this chapter first we give a brief introduction to 
linear sequential circuit (LSC) in general. In Section 4.2 we 
discuss under what conditions the response of a composite sequ- 
ence over GF(2) or any of its extensions can be obtained by the 
sum total of responses due to constituent sequences. In Section 

4.3 we discuss a method of obtaining a set of basis functions 
for the vector space of periodic Galois functions characterized 
by difference equations. This method is made vise of in Section 

4.4 for obtaining basis functions for the vector space of 
periodic Galois sequences of period equal to the cycle length 
of the scrambler, which when given as input to the scrambler 
results in an output sequence of the same period, irrespective of 
initial state. Xn Section 4.5 we give a method to obtain 
Fourier series expansion of periodic binary sequences over 
extension fields of GF(2) . In Section 4.6, we discuss concept 

of eigensequences and eigenvalues of LSCs . This concept is made 

use of in Section 4.7 to carry out spectral analysis at input 

and at output of 2,3, 4 and 5 stage scramblers for various 

cases of input and output periods, taking eigensequences appro- 

2 3 4 5 

priately from extension fields GF{2 ), GF(2 ), GF(2 ) and GF(2 ). 
Some general observations are given at the end of the section . 
Operation tables for various extension fields are given in 
Appendix a. 
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4.1 Linear Sequential Circuits (LSC) 

A block diagram of a LSC is shewn in Figure 4.1. 




Essentially it is made up of delay elements which serve as the 
memory and logic networks constructed from logic elements which 
perform the Galois field operations of addition and multiplication 
by a constant. An input sequence { u^,| applied to a LSC results 
in an output sequence | y-^ which is a linear function of the 
present input value and the present state. The present state 
in turn is a linear f -unction of past states and past inputs. 

A linear continuous system is characterized by a set of 
first-order differential equations. In the case of a LSC, it 
is characterized by a set of first-order difference equations. 


These equations can be written in matrix form as 

Xjc+1 - *** + Bu k (4ai 

Vic = cx k + ^ (4 - 11 

where, X is state vector, u and y are input and output vectors 
respectively. For single input single output IEC u and y are 
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both one— tuples. The characterizing matrices A, B, C and D are 
finite field valued and of consistent dimensions . The matrix; A 
governs the autonomous behaviour of the LSC and is called 
characteristic matrix. 

Let X Q denote the initial condition or state of the 
machine at k = 0. Then, it can be shown that 

c “X + 



CA k ~j -1 


BU. 




(4.2) 


We can see that the response of the system has two parts, one is 
zero input response and the other is zero state response. Repre- 
senting these by £ y z1 and |y Zs respectively, we can write 


y Zi,k + y Zs,k 


where 


y Zi,k = A ^o 

k-1, . 1 

and y Zs,k = ^ CA BUj + LUj, . 


(4.3) 

(4.4) 

(4.5) 


4.2 Total Response from Responses to Component Sequences 

Let the input sequence | be expressed as a bit by 
bit sum of its component sequences, i.e. 




- 2^ u ± i 

i=i 1/] 


then let us investigate under what conditions the output sequence 
is given by 

r, 

y k = 2_, y i/k 
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where the sequence j^y^ ^ is the output for the input sequence 
£ u r alone. Let 


X l,0 

X 2,( 

Then from Equation 

x 14 

= 

X 1 , 2 

= a 2 x 


r,0 


* 1,0 

*.Bu, A + Bu. . 

1 f 0 1/1 


x i,k - A \ * A *~X,o + 


+ ABu l,k-2 + Bu l,k-1 


Similarly, 

X 2,k 


~ A^X _ 4- A Bu« A + 


. k-1 

-o A bu 2,0 


+ ABu 2,k-2 + Bu 2,k-1 


X . 
r,k 

By adding 


“ A \ + ^'X.O + 


+ ABu r,k-2 + Bu r,k-1 


xr~i 


s X. . = A^X +. ..+ A*X_ + ' k_1 

m l '* ° 


A - Bu o +. . .+ ABu k _ 2 + BU ]< _ 1 


(4.6) 


r times 

Also, for the input sequence we can write 

*k - + aK ~ 1bu o * ■■■ + ABU k-2 + Bu k-1 (4 - 7) 

Prom Equation (4.6) and Equation (4.7) it follows that 
r , 

SX* = ** i£ (1) x ° “ - or . 

(ii) rA^X = A*X o 


(4.8) 
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The second condition is satisfied for any odd number r since 
the operation is of characteristic 2? i.e. if the input is 
split into odd number of component sequences. Now for each 
tJle output can be written as. 


■ i,k 


CX + du. . 

i,k 


r. 


f?r 


c x i,x + D 


c 2_. x i k + Dir 


(4.9) 


Also for input sequence f , we have 


f K = °Sc + Du k 


(4.10) 


From Equation (4.8) and Equation 


*k 


r 


f=T i/k 


(4.10) we can say that 


if = 0 or the input sequence is split into odd number of 
component sequences otherwise to the summation of individual 
responses, we must add the aero input response to get. the total 
response y,., i.e. 



if r is even. 


(4.11) 


4.3 Basis Functions for the Vector Space of Periodic 


Galois Functions J_ 2 


In a continuous linear system characterised by differ- 
ential equations, various signals in the system are represented 
by losing functions such as e , e sinut and t e cos^j t. 


There is a corresponding set of functions which are used to 
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describe the signals in linear machines . One such function is 
the exponent function as defined below. 

Let (P be any non -zero element of GF(p n ) . Then the 

function 

e(t ) = p 1 

which is defined for all integral values of t, is called the 
exponent function. If p is of order k, then 

p t+lc = P 

Hence e(t) has period k. For example if p = a 5 where a 5 € GF(2 4 ) 
then 

e{t) = (a ) = 1 a a la a ... for t = 0,1/... 

The response of a system to an exponent function is called its 
exponential response. 

Periodic Galois Functions: Any sequence gjt) defined for 

integer values of t and taking values from a Galois field 

GF(p ) is called a Galois function. The linear delay operation 
apply to these functions. The delay operation is indicated as 
D^g ( t ) = g(t-i) where D is referred to as the delay or D operator. 
If g(t ) is such that g(t) = g(t-T), then the smallest T for 
which this is so is called the period of g (t ) . 

Periodic Galois functions are not amenable to represen- 
tation in a simple closed form. Instead, these functions can be 
defined in terms of the solution to linear difference equations . 
The general form of such an equation can be written as 


g(t) + a^t-l) + 


+ a q g(t-q) 


0 


(4.12) 
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where a^' s are selected from GF(p n ). Using the delay operator D, 
this equation can be written as, 

[1 + a x D + ... + a q D q ] g(t ) = F(D)g(t) = 0 (4.13) 

The polynomial F (D) is referred to as a delay polynomial, 

and a function g(t), which identically satisfies the equation 

above for all values of t, is said to be a solution to the delay 

equation. For example any periodic function with period T will 

T 

satisfy the equation [1 - D ] g(t) = 0. 

The set S ? of all such g(t) which satisfy Equation (4.13) 
form a q dimensional vector space over GF(p n ). This space is 
invariant under linear delay operator D. Any function g (t ) in 
this space can be written as 


g(t) 



b ± *i<*> 


(4.14) 


where ^ (t ) , 0 2 <t), 0^(t ) form a set of q linearly indep- 

endent functions; i.e. they form a basis for Sp. Now let us 
see how these basis functions can be obtained. 

Any polynomial F(D) can be written as the product of 
its prime factors, i.e. 

F(D) = [^(D)] 1 [p. 2 < D) 3 2 ***fP r (D) 3 r (4.15) 


where \i. (D) is an irreducible polynomial of degree q^. Each 
'© . 

factor [ ii^(D)] 1 corresponds to a linear subspace and associated 
with each, is a set of q^e.^ basis functions 



j = 1 , 2,...,.q ± / k » 1,2, 
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spanning a vector space s{V e^. If 1 (k) 0 . .(t)l and 

((k) 1 ^*'i' '■* ^ x, 3 J 

f 0 g, j (t) V are the functions corresponding to factors 

and [n (F)] g respectively, then the basis associated 
with F(D) = [^(D)] i [(i g (D)] eg is [ (k) 0 io .(t)^V £ (u) ^ g ^(t)]. 
Thus we can find a way to obtain basis for each subspace 

S [ji.(D)] e i t ^ ien ‘union’ of all such basis gives a basis for the 
vector space S p . 

0 » 

To find the basis functions associated with [^(DjJ 1 
first select q^ factions that satisfy the delay equation 
p.^(D) g(t) =0. Next q^ functions are selected so that they 
satisfy 

[^(D)] 2 g(t ) = 0 , but not ^(D) g(t) = 0. 


i, T_ 

This is done successively to obtain e i set of q^ basis 

functions so that they satisfy 

0 0 
[H ± (D)] 1 g(t) = 0 but not [p ± (D)] 1 " g(t) = 0 


Let T q be any 


Now consider [n-(D)] e = 1 + a.D + ... + a D q 

q T 

value of T such that [ p,(D)] e divides [l — D °] * Then the period 

of every element in S will have a period which is a 

[li(D)] S 

divisor of T . To generate the basis function of S- • 

° MD)] e 

first generate the function h(t) such that 


h(t ) 


1 

0 


t = nT Q n = 0,1,... 
otherwise 


Oi 

which satisfies the equation [l - D ]h(t) = 0. 


(4.16) 

Using h<t) 


( 1 ) 


0 1 (t ) can be generated as 


the function 
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( 1 ) 


0 x (t) 


(1 — D °) 
\ i(D) 


h(t ) 


(4.17) 


Note that the f motion is not zero and that y(D) ^ V-^t ) = 0. 
Hence it is selected as first basis f motion. If q = 1 this 
will be the only basis function. Otherwise other q-1 basis 
functions are given by 


<1) 0 i (t) = D 1 ” 1 [ (1) 0 1 (t)] , i = 2 q (4.18) 


In a similar manner the q basis functions associated with 
[ |i(D)] ^ for j = 2,3,...,e are given by 


<J> 


^(t) 


= D' 


i-1 


T 1 


1 - D 


h(t) 


i = 1,2, ... ,q (4.19) 


After obtaining basis functions associated with each prime 
factor, the basis for the vector space is obtained as the 
union of all these basis functions. 


Example 4.1: 

Consider a function 

g(t) such 

that 


F(D) g(t ) 

2 3 

= (1 + + D J + 

D 4 ) g(t) 

= 0 

We can write F(D) as 




F(D) = 

(1 + D) (1 + D + D 3 ) 

* 


Basis 

Associated with 1+D: Since 

(1 + D) is a divisor 

(1 + 

D) itself 





h(t) = 

11111111... 



and 

{1> 0 1 (t) 

1+D -u i*. \ _ 

- m h(t) ■ 

h(t ) 


i.e. 

0 x (t) = 

11111111... 
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Basis Associated with 1 + D + d 3 : The factor (1 + D + D 3 ) 

7 

divides 1 + d 

h(t ) - 1000000100..., 

7 

and = — -- - - - — =■ h(t ) = (1 + D + D 2 + D 4 ) h(t) 

1 + D + 

= h(t ) + h(t-l) + h(t-2) +h(t-4) 
i.e. ^^(t) = 1110100... 

7 

{1) j2L(t) = d — L±,R — _ h(t) 

1 + D + D 

= (D + D 2 + D 3 + D 5 ) h(t) 

= 0111010... 

and = D 2 — - — — 3 h(t) 

J 1 + D + D 

= (D 2 + D 3 + D 4 + D 6 ) h(t ) 

(1> 0 3 (t) = 0011101... 

Hence we have obtained a basis for the space and the basis 
functions are, 

0^(0 = 1111111... 

0 2 (t) = 1110100... 

0 3 (t) = 0111010... 

0 4 (t) = 


0011101 ... 
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4.4 Application to Scramblers 

New we apply the method of the previous section to 
three and four stage scramblers for obtaining input sequences of 
periodicity equal to cycle length q and giving rise to output 

sequences of the same period. 

# 

Illustration 4.1: If Figure 4.2, a three stage scrambler with. 

3 

characteristic polynomial X + X + 1 = 0 with a periodic input 
of period 7 is represented. 



Output y(t) 


3 

Figure 4.2 Three Stage Scrambler: c(x) = x + x + 1. 


rp 

If the initial state vector is Z = [u^u^ . . . u^X^X^^] , then 
the state vector at any instant X is given by 


where T is the transition matrix given by 

0000001000 
1000000000 
0100000000 
0010000000 
oooioooooo 
r ooooiooooo 

0000010000 
0000001011 
0000000100 
000000001 0 
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7 

Prom this , T can be found, to be given by 

1000000000 

0100000000 

0010000000 

0001000000 

7 0000100000 
T = 

0000010000 

0000001000 

1011100100 

0101110010 

0010111001 


1 2 

Sj- ee the output y(t) is given by, y(t) = [ T Z] g, [T z] g, .../ 

[T z]g ... where [Z] n stands for value at n location of the 
vector Z. Hence it follows that for the output to be of period 7, 


7 7 

Z = T Z = Z . 

If this condition is to be satisfied irrespective of initial 
state of the scrambler/ the input sequence u(t) must satisfy 
following equations/ i.e. 

u^ + Uj- + + u 3 = 0 

u, + u ' + u 0 + u~ = 0 (4.20) 

6 4 3 2 

u 5 + u 3 + u 2 + “l = ° ' 

Writing in the delay polynomial form Equation (4.20) becomes. 
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(1 + D 2 + D 3 + D 4 ) u(t) = 0 

D(1 4- D 2 + D 3 + D 4 )u(t) = 0 (4.21) 

D 2 (l + D 2 + D 3 + D 4 )u{t) = 0 

Any solution to above equations is an input sequence u(t) of 
period 7 resulting in output period of 7 itself. As discussed 
in Section 4.3, all such input functions u(t), form a vector 
space S F ( D ) where 

F(D) = 1 + D 2 + D 3 + D 4 (4.22) 

We have obtained a basis for this space in Example 4.1. The 
basis functions are 

0 x (t) = 1111111... 

JZL(t) = 1110100... 

(4.23) 

0 3 (t) = 0111010... 

0 4 (t) = 0011101... 

Note that the basis functions 0 3 (t) and 0 4 (t) spans the 

3 

vector space of 2 sequences generated by the 3 stage M— sequence 

3 2 

generator with characteristic polynomial c(x) = x +x +1. 

4 

Also it can be seen that (a) there are (2 -2) =14 sequences 

excluding all zero and all one sequences of period 7 giving rise 

to output period of 7 itself irrespective of initial state of 

the scrambler, (b) these 14 sequences are M-sequences generated 

3 2 

by the structure with characteristic polynomial c(x)=x +x+l 
and their complements . These facts confirm observations made 


in Section 3.7. 
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Similarly for a scrambler structure of characteristic 
polynomial c(x) = x + x^ + 1, it can be shown that any input 
sequence u(t) of period 7 and giving rise to output period 7 
itself irrespective of initial state, must satisfy the equation 

F(D> u(t) = (1 + D) (1 + D 2 + D 3 ) y (t ) = 0 

A basis for the vector space formed by such functions can be 
shown to be 

^(t) = lllllll... 

0 2 (t) = 1011100... 

(4.24) 

0 3 (t) * 0101110... 

0 4 (t) = 0010111... 

M 0+-0 that the basis functions , 0 3 (t) and 0 4 <t) spans the 

3 

vector space of 2 sequences generated by the 3 stage M~sequence 

3 

generator of characteristic polynomial c(x) = x + x + 1. Also 
it can be seen that 

(a) there are 14 sequences excluding all zero and all one 
sequences of period 7 giving rise to output period 7 
itself, irrespective of initial state of the scrambler, 

(b) these sequences are M-sequences generated by the structure*, 

3 

of characteristic polynomial c(x)=x +x + l and their 
complements, which is in confirmity with the observation 
made in Section 3.7. 

Illustration 4.2: Consider a 4 stage scrambler of character- 

istic equation x 4 + x + 1 = 0 and with a periodic input of 
period 15 as depicted in Figure 4.3. 
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Output 

Figure 4.3 Four Stage Scrambler; c(x) = + x + 1. 


Proceeding exactly like in Illustration 4.1, we can 
obtain following equations to be satisfied by an input sequence 
u(t) to give rise to output sequence of period 15 itself irres- 
pective of initial state of the scrambler 

U 1S + u 12 + u U + u 9 + "7 + + u 5 + u 4 = 0 

u 14 + U 11 + u 10 + u 8 + u 6 + u 5 + u 4 + u 3 “ 0 

(4.25) 

u 13 + u 1Q + u g + xij + u 5 + u 4 + u 3 + u 2 = 0 

u 12 + Ug + u g + u 6 + u 4 + u 3 + u 2 + Uj - 0 


Writing this using delay operator, we get 

(1 + D 3 + D 6 + D 8 + D 9 + D 10 + D 11 ) u(t ) = 0 

D(1 -f D 3 + D 6 -f D 8 + D 9 + D 10 + D 11 ) u(t) =0 

(4.26, 

D 2 (l + D 3 + D 6 + D 8 + D 9 + D 10 + D 11 )u(t) = 0 
D 3 (l + D 3 + D 6 + D 8 + D 9 + D 10 + D 1X )u(t) = 0 

As discussed in Section 4.3, all such input sequences u(t) form 
a vector space S p (D), where 
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F(D) - 1 + D 3 + D 4 + D 6 + D 8 + D 9 + D 10 + D 11 

F(D) = {1 + D)(l + D + D 2 )(D 4 + D 3 + D 2 + D + 1)(D 4 + D + 1) 

(4.27) 

Now we will find basis functions associated with each factor. 


Basis 

for (1 


D) : Prom Example 4.1, we have 


0 x (t) 


11111111... 

Basis 

for (1 


D + D 2 ): (1 + D + D 2 ) divides (1 + D 3 ) 

# * 

h(t ) 

S 

1001001001... 

*» 

0 2 (t) 

S 

—= — — 5 Mt) 

1 + D + D 



SK 

(1 + D) h(t) = 110110110... 

and 

0 3 (t) 


D0 2 (t) = 011011011... 

Basis 

for (1 

4* 

D + D 4 ): ( 1 + D + d 4 ) divides 1 + D 1 5 . Hence 


h(t ) 

=5 

100000000000000,100. . . 

♦ 

* • 

**4 (t) 


i _15 

h(t) 

1 + D + D 



5= 

1 + D + D 2 + d 3 + D 5 + D 7 + D 8 + D 11 h(t) 


0 4 <t) 

5= 

111101011001000. .. 


%<t) 

= 

D0 4 (t) = 011110101100100... 


0 6 <t) 

S 

D 2 0 4 (t ) = 001111010110010... 


0 7 (t) 

S 

D 3 0 4 (t) * 000111101011001... 


000111101011001 ... 
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Basis for <D 4 + D 3 + D 2 4 - D + 1): (D 4 4 D 3 + D 2 + D + 1) 

5 

divides 1 + d 


h(t) = 10000 , 100 . . . 

Basis function associated with this are. 


0 8 (t) = --4 h(t) 

D 4- D 4* D 4-04-1 

= (1 + D) h(t ) 


0g(t) 

J2 9 (t) 

^io (t) 

0n<t> 


= 1100011000 ... 

= T& g(t) = 01100,01100... 

= D 2 0 g (t) = 00110,00110... 

* D 3 0 8 (t)= 00011 , 00011 ,... 


Hence, a set of basis functions for s F ^ D j are, 

0 x (t) , 0 2 (t), ..., 0 lo (t ) and 0 n (t). 

Prom this, it can be seen that 

(a) there are 2 ^ - 2 3 - 2 + 2 “* - 2 -2 = 2010 sequences 

of period IS which gives rise to output period of 15 
irrespective of initial state, 

(b) there is a subspace spanned by functions 0 1 (t), 0 4 (t), 
0 5 (t), 0 6 <t) and 0rj (t ) consisting of all sequences gener- 
ated by the M-sequence generator with characteristic 
polynomial c(x) = x 4 + x 3 + 1 , and sequences complement 
to them, which is in confirmity with the observation made 
in Section 3.7. 
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Similarly/ for the scrambler with characteristic polynomial 
4 3 

c(x) = x + x + 1, we can get a condition that is to be satis- 
fied by an input sequence of period 15, to give rise to output 
period 15 irrespective of initial state. The condition is, 

F(D) u(t) = (1 + D)(l + D + D 1 2 )(D 4 * + D 3 + D 2 + D + 1) 

(D 4 + D 3 + 1) u(t) = 0 (4.25) 

We have already obtained basis function associated with first 
three factors earlier. Hence let us obtain basis functions 

A -5 4 "a I c. 

associated with (D + D +1). (D + D +1) divides 1 + D - 
Hence h(t) can be taken as, 

h ( t ) = 100000000000000 ,100... 

and a set of basis functions can be obtained as 


0 4 (t) = 


1 -f D 


15 


_4 , _3 , . 
D + D +1 


h(t ) 


1 + D 3 + D 4 + D 6 * + D 8 + D 9 + D 10 + D 11 h{t) 


0 4 (t) = 100110101111000... 

0 5 ( t) = D0 4 (t) = 010011010111100... 

0 (t ) = D 2 0.(t) = 001001101011110... 

0 ? (t) = D 3 0 4 (t) = 000100110101111... 

These four functions and functions 0 1 (t) , 0 2 (t) , 0 3 (t) , 0 Q (t) , 

0 g (t ) , 0 lO (t) and ^n (t) obtained earlier constitute a basis 

for the vector space of all such functions u(t ) satisfying 

Equation (4.28). From this it can be seen that 
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(a) there are 2 11 - [ 2 3 - 2 + 2 5 - 2] - 2 = 2010 sequences of 
periodicity 15 resulting in output period 15 itself irres- 
pective of initial state, 

(b) there is a subspace spanned by basis functions 0^(t), 

0 4 (t), 0 5 <t), 0 & (t) and 0 7 (t) consisting of all sequences 
generated by the M-sequence generater with characteristic 
polynomial c(x) - x^ + x + 1, and sequences complement to 
them, which is in confirmity with the observations made in 
Section 3.7. 

4*5 Fourier Series Expansion of Periodic Binary Secruences 
Over Finite Fields ” 

A periodic real or complex signal can be expressed as 
an infinite sum of exponential functions in complex field. 
Similarly in the case of a periodic function over GF(p) it is 
possible to express this as a finite sum of exponent function 
over an appropriate extension field GF ( p n ) . The difficulty in 
this case is that it is possible to do so for periodio. ^signals with 
certain periodicities only. Only such functions having a 
periodicity i, where i is a divisor of p 1 -! can be expressed 
in this form. For sequences over GF(2) satisfying above const- 
raint, we give below a method for obtaining Fourier series 
expansion. 

p m -l 

From Equation (2.5) it is clear that (x - 1) splits 

into (p - 1) linear factors in GF(p ),i.e. 

m mi 9 

(x p “ 1 - 1) = (x - a°)(x - a) ... (x - a p ) 

Hence all non-zero field elements are various roots of unity. 

Now consider a matrix M whose columns consist of exponent 
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« rs ll ' «'*) 

functions of non-zero elements a to a;P ~ . The matrix M has 
following properties. 

(a) It is a square matrix of dimension (p m -l) x(p m -l) 

(b) H is a symmetric matrix since 

m ii ” = (a^)^ = m.,, i = 0,...,p m -2 

ji 

j = 0,...,p m -2 

(4.29) 

(c) For each column j , there is a column k such that 

m ik = m ij because, for every element a-* there is an 
v 

inverse a and if that is so 
fV = [(a 3 )" 1 ] 1 = [(a 3 ) 1 ]" 1 


JL « 0 + 


m ik = m ij 


(4.30) 


(Note that the first column is all ones and hence it is its own 
inverse). Also for such columns j and k 


m. 


1=0 


m . 
p -2 




ik' '“‘ij 


i=a 


mod p. 


(= 1 in GF(2 m ) ) 


(4.31) 


(d) Multiplication of respective elements of any two columns 
give another column, i.e. 

(m. . ) (m. v ) = (ah 1 (a*) 1 = (a^ 5 ^) 1 = (a 1 ) 1 = m 

(4.32) 

and for any two columns j and k not satisfying the relation 
in (c) 
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p m ~2 
i=0 1 ' 




0 (4.33) 


From the above properties, it follows that the matrix 
can be obtained by interchanging appropriate columns of M. 


Illustration 4.3: The matrix M given below is a matrix having 

exponent function in GF(2 3 ) as its columns. 



It can be seen that the pairs of columns (2, 7), (3, 6), (4, 5) 
satisfy property (c). Hence can be written as. 



Now let us say a periodic binary sequence of period i (such 
that ±L is a divisor of 2^-1 ) can be expressed as 



(4.34) 
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where £ a J denotes exponent function of element cf^. This can 
be written in matrix form as 

M 2 = A (4.35) 

where A - [a^ a m ] T , C = tc^,...^ m 3 T 

2 -2 2 “2 

and the coefficients c Q , . .., c can be obtained as 

2 “ 2 

C = M A (4.36J 

Once we are able to find the Fourier coefficients, the 
next question which arises is that whether there can be any 
well defined, exponential response for a given LSC due to each 
one of these exponent functions. Because, in that case the 
response due to any binary sequence (with appropriate perio- 
dicity) can be obtained from exponential responses due to its 
constituent exponent functions. This aspect is discussed in 
the next section. 

Example 4.2: let us now find out the Fourier coefficients 

(spectrum) of the periodic sequence £ u^J = 1001011,1001011... 

Let the coefficient of any [ e^J be denoted by U^, 

3 

where e 6 GF(2 ). Then these coefficients are obtained as 


where , 

Se 

■ £ U 1 U a V ••• V 1 

and 


- tv U 1 ••• U 6^ T 
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i .e . 


or 


t? 


1 

1 

i 

1 

1 

1 

1 


V 

1 

a 6 

a 5 

a 4 

o 3 

a 2 

a 


0 

1 

a 5 

a 3 

a 

a 6 

a 4 

a 2 


0 

1 

a 4 

a 

a 5 

a 2 

6 

a 

a 3 


1 

1 

a 3 

a 6 

2 

a 

5 

a 

a 

4 

a 


! 0 

1 

a 2 

a 4 

6 

a 

a 

a 3 

5 

a 


1 

1 

a 

a 2 

a 3 

a 4 

a 5 

a 6 


1 

w 










U^_ = l+l+l+l = 0 

4 2 

U a = 1+a + a + a = 1 

U 9 = 1 + a + a 4 + a 2 = 1 

or 

U » 1 + a 5 + a 6 + a 3 = 0 

or 

U = 1 + a 2 + a + a 4 = 1 

or 

U = l + a 6 + a 3 + a 5 = 0 

or 

U . = 1 + a 3 + a 5 + a 6 = Q 


4.6 Biqensiqnals and Eigenvalues of LSC 

We will now study ‘Spectral theory* on finite fields 
for LSCs as suggested by Sangani [13]. 

For an operator M on a vector space, if there exist a 
scalar \ and a non-zero vector u such that 

MU = >.U (4.37) 

then X is said to be an eigenvalue and U, the corresponding 
eigenvector. This equality which is very familiar with real 
or complex matrices on a finite dimensional space also holds 
for an arbitrary operator on an abstract vector space. This 
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fact is made use of to define the eigenvalues and eigenvectors 
or eigenfunctions of a single input -single output LSC by 
viewing it as a linear operator on the infinite dimensional 
vector space of ’two-sided* sequences with values in GF(p) . 
Hence Equation (4.1) can be written as 


*k+l 

= **k 

+ 



xu k 


+ Du k - 

- CD < k < CD 

(4.38) 


where { u^ is an eigenfunction and A the eigenvalue,, both taking 
their values in some extension field of GF(p). Any extension 
field of a field plays essentially the same role for an LSC 
defined, over that field as complex numbers do for an analog 
circuit defined over the real field. 

Consider a periodic exponent function £ e' K j", k = 0 , +1 , 
+2.../ as input to an LSC. Where e is an element in some exten- 
sion field of GF(p) over which LSC is defined. Transfer function 
S(e) of this LSC is defined as 

S(e) = D + C(el - A)" 1 B (4.39) 

Here S(e) is a rational function in e with coefficients 
in a finite field, whose poles and zeros are well defined elem- 
ents of the extension field. For each e, except for the cases 
where (el - A)” 1 is not defined (does not exist), the sequence 
“ od < k < <od, is an eigenfunction of the LSC with 
eigenvalue S(e). The eigenvalues of LSC are elements of the 
extension field and the eigenfunctions take their values in the 
same extension field, just as the eigenvalues and eigenfunctions 
defined over the real field are complex. 





The following theorem gives the condition under which 
S(e), the transfer function defined earlier and the function 
{e i are respectively an eigenvalue and corresponding eigen- 
sequence of the LSC defined by Equation (4.38). 

Theorem: Let 

Xfc+i * ^ + Bu k (4.40) 

Y k - CXj, + DUj, , - cd < k < CD (4.41) 

characterize a LSC over a finite field. Then for each e in an 
extension field for which 

S(e) A c(le - A) -1 b + D 

is defined, there exists a unique state X q _ X(e) such that the 
sequence 

= { e k | , k = 0, +1, +2, ... 

is an eigensequence for the LSC defined by Equation (4.40) and 
Equation (4.41). The corresponding eigenvalue is S(e). The 
unique state X Q is given by 

X Q = X(e) = (Ie - A)” 1 B . (4.42) 

It can be shown that 

Xj^. = (Ie - A) 1 Be k (4.43) 

is the 'unique solution to Equation (4.40) with X q given by 
Equation (4.42). Substituting Equation (4.43) in Equation (4.41), 


we get 
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y k * CX k + 

* [C(Ie - A)~ a b + D] e 
i.e. y k = [S(e)]e k 


k 


<4.44) 


(4.45) 


Let 'us now consider an input sequence .^u^ which is a 
shifted version of an eigensequence [e k j; i.e. 

f ^ j = b h k i ' b 6 GF(2 n ), b/0 

Then from Equation (4.43), 

\ = b(Ie - A)" 1 b e k 

Substituting this in Equation (4.41), 


(4.46) 


Y k = b [ C(Ie “A)” 1 B + D] e k 

= S(e) b e k = S{e) (4.47) 

i.e., the coefficient of {e k ] at the output is eigenvalue times 
the coefficient at input . We can observe that the function ^ u^j. 
is an eigensequence with S(e) as eigenvalue and the unique state 
is given by 


X(e, b) * b(Ie - A)” 1 B = bX(e) 


(4.48) 


where X(e) is the unique state X Q such that the sequence e 
and S(e) are respectively an eigensequence and corresponding 
eigenvalue of the LSC. 

Now consider a sequence { u^ given by Equation (4.45), 
given as input to a LSC operating in GF(2) and whose initial 
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state is 2 $ro # The total response ^ ^ can be written 


as 


_£V — V * 1 

% = ZEI ^ ~ :i ” Bu jt + Dv^ + A' k x(e / b} + A^XCe^b) 

From Equation (4.47 ) this can be written as 

y k = S(e) u^ + A^Ce.b) (4.49) 

If initial state is non-zero, say X , then 

a 

y k * S(e) u k + A k X(e / b) + A^x (4.50) 

where first two terms give zero state response and the third 
term gives' zero input response. 


4.7 Application of Spectral Theory to Scramblers 

In this section we will apply spectral theory to 2, 3, 

, k r 

4 and 5 stage scramblers . Eigensequences £ e $ are taken from 
extension fields. GF (2^ ) , GF(2 3 ) , GF(2^) and GF(2^) in each case. 
Spectral analysis at input and at output of the scrambler for 
some class of periodic input sequences is carried out and 
observations made thereof are summarised at the end of the 
section. Throughout this section the coefficient of an eigen- 
sequence je k | at input and output are denoted by U g and 
respectively. 

4.7.1 Two Stage Scrambler: 

The two stage scrambler with characteristic polynomial 
c(x) = X 2 + X + 1 is shown in Figure 4.4. 
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Figure 4.4 Two Stage Scrambler; c(x) = + x +1. 

The characterising equations can be written as 

\+l = AX k + Bu k 

(4.51) 

Y k = 

where 

1 1 1 1 j 

A = , B = , C = [l l] and D = 1 

1 0 I L° . 

2 3 4 

Let us new consider eigen sequences from GF(2 ), GF(2 )<» GF(2 ) and 
5 

GF(2 ) and we will find out spectral coefficients corresponding 
to each of them at the input and at the output. 

2 

(a) Eigensequences from GF(2 ); For e = 1, the Het(Ie - A) is 
given by 

0 1 I 

Det(Ie - A) = =1 

1 1 

Therefore for e — 1 eigenvalue is defined and the unique state 
X Q = X(l) is given by 

X(1 ) = (le - A)” 1 B = 

S(l) = C(le - A)" 1 B + D = 1+1+1 = 






113 


For e = cc. 


Detde - a) 


1+a 1 
1 a 


Therefore eigenvalue is not defined 
sequence of this LSC. 



is not an eigen- 


For e = a 


2 




Det(Ie - A) 


1 + c? 
1 


1 

a 


= 0 . 


Therefore eigenvalue is not defined 
sequence of this LSC. 

Now if the sequence 111... is given as input with the 
state X(l) then, output sequence is also 111... i.e. the Equation 
(4.44) is satisfied. For two stage scrambler there is only one 
case wherein the output is of period 3 irrespective of initial 
state,* Spectral coefficients at input and at output for this 
case and different initial conditions are obtained as explained 
in Section 4.5 and are listed in Table 4.1. Spectral coefficients 
for zero input responses are also given in the table. 

Note that * 

(i) the spectral coefficients of eigen sequences and 

a ^ j for which eigenvalue is not defined, are zero at 
the input. 

(ii) For e = 1 where eigenvalue is defined, 

Y x = S(l) U x 


and 


i 2k? . 

r r 


is not an eigen- 


* State at instant k = 0 
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Table 4.1 

Spectral Coefficients 
Sequences of Period 3? 

at Input 
2 -stage 

and Output for 
Scrambler 

Initial 

States 

X 1 X 2 

Input/ Output Sequences 

i— — - ~ 

Input/Output Spectral 
Coefficients 

vv w VV 



in,... 

l 

0 

0 

0 

0 







100,... 

l 

1 

1 



111,.** 

l 

0 

0 

0 

1 







010,... 

l 

a 2 

a 



111,... 

l 

0 

0 

1 

0 







001,... 

l 

a 

2 

a 



111,... 

l 

0 

0 

1 

1 







111,... 

l 

0 

0 

0 

1 

■ - 


- 




110,... 

0 

a 

a 

1 

0 

- 


o 




101,... 

0 

£* 

a 

a 

1 

1 

- 


- 




Oil,... 

0 

i 

1 
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i) e new coefficients which appear in those places where 

not 

eigenvalues are/defined, are due to the sequence generated 
by the state given by 

*T * U 1 X(1 > + x a (4.52) 

where X & is actual initial state of the scrambler. This 
is justified by Equation (4.50) . 

Example 4.3: In the case when the input is 111... and initial 

state is 01, X T is given by 



and the sequence generated by this is 101,101.... The spectral 

coefficients of this are Y. = 0, Y = a 2 and Y « = a. From the 

la q2 

Table 4.1 it can be seen that these are the new coefficients 
which appear at the output. 

(b) Eigensequences fromGF(2 3 ): For the 2-stage scrambler, 

eigenvalues S(e), for each eigenfunction £ e k \ and their unique 
states X(e) are given in Table 4.2. In this case eigenvalues 
are defined for all eigensequences. 

Now let us consider 7 basis functions which span, the 
vector space of all 7-tuples as input. Each one of them gives 
rise to an output sequence of period 7 when the scramble is in 
critical state. The spectral coefficients at input and output 
for s uc h cases are listed in Table 4.3. It can be seen that 
(i) For any eigensequence {e k |, its coefficient Y Q at output 
is given by 
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Table 4.2 


Eigenvalues and Unique States in GF(2 3 ); 2 -stage 

Scrambler 




where X„ is the critical state for that sequence, 
c 

This can be justified as follows. From Equation (4.11) 
and Equation (4.50) the output response is given by 


U S(e)e k + A k / , U X(e) + A^X (4.54) 

e 

In this, the response due to second and third terms is 
of period 3. Hence for the output period to be 7, this 
response should be all zeros. This will be so only if 
condition (ii) is satisfied- 
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Table 4.3 Spectral Coefficients at Input and Output for 
Sequences of Period 7; 2-stage Scrambler 


* 

Critical « 

States X ! 

c , 

Input/Output Sequences 


Input /Output 
Coefficients 

Spectral 

X 1 

x 2 ; 
— .... t. 









1 

i 

1000000,... 

1 

1 , 

1 

i 

1 

1 

1 



1011011,... 

1 

a 4 

a 

a 

J2 

cc 

a 4 

a 2 



0100000 








1 

0 

1101101,... 

1 

a 3 

a 6 

a? 

C£ 5 

£ 

d 3 



0010000,... 








0 

1 

mono,... 

1 

i 

a 4 

a 2 

a 

a 

a 4 



0001000, . . , 








1 

1 

onion,... 

1 

a 

a 2 

a 6 

a 4 

cc 3 

a 5 



0000100,... 








1 

0 

1011101,... 

1 

1 

1 

a 3 

1 

a 5 

i 

1 

t 



0000010, . . . 








0 

1 

1101110,... 

1 

a 6 

5 

a 

1 

„3 

CC 

1 

1 

— -ry - 



— _ 



— 



— 



- — — 


0000001 , 


1 a 5 a 3 a 4 a 6 a 2 


1 


0 


0110111 ,... 


From Section 3,9 it follows that any input sequence of 
period 7 , when the scrambler in critical state, will give rise to 
an output sequence which is seme linear combination of the 7 
output sequences given in Table 4.3. The c or re s ponding critical 
state is given by the same linear combination of critical states 

listed against each output sequence. Therefore the above two 

•> 

properties hold good in all cases where input and output periods 

\ 

are 7 . 


Example 4.3: From Table 4.3, for the sequence 0100000,01. 

we have 


\ 1 * 

\ U X = 
2j I e e 


1 

6 

a 3 

5 

~* 6 ~ 

4 


m 

+ a 

2 

+ a 

4 

+ a 

2 . 

1 


l CC 


a 


foe 

L i 


a 5 ' 

2 

R ' 
O' 

J 


a 3 


+ a 


+ a 




a j 


a 4 



i 


" 1 ’ 


•, 3 

1 + a + a 


0 


- 




the critical state for that sequence. 


* / 


4 5 

(c) Eigensequences from GF(2 ): In this case except for e = a 

and a 10 for all other eigensequences, eigenvalues are defined. 

Eigenvalues for all eigensequences and their unique states are 

listed in Table 4.4. We will do the analysis for the following 

cases. 

(i.) Input period is 5 and output period is also 5 
<ii) Input period is 5 and output period is 15 
(iii) Input period is 15 and output period is 15. 
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Table 4.4 Eigenvalues and Unique States in GF(2 4 )? 2-stage 
Scrambler 


e 

T 

1 X(e ) 

j X, 

- 1 - . 

S 

« 

1 

X 2 j 

S(e) 

i 

1 

1 

1 

a 

a 6 

a 5 

a 7 

a 2 

a 12 

a 10 

a 14 

a 3 

a* 0 

a 7 

a 13 

a 4 

a 9 

a 5 

a 13 

a 6 

a 5 

a 14 

a 11 

a 7 

a 3 

a 11 

a 10 

a* 

a 3 

a 10 

a 31 

a 9 

a 5 

a 11 

a 14 

a 11 

a 9 

a 13 

a 5 

12 

io 

13 

7 

a 

a 

a 

a 

a 13 

a! 12 

a 14 

10 

a 

14 

6 

7 

5 

a 

a 

a 

a 


In Table 4.5 we have listed spectral coefficients at input and 

at output for 5 basis functions of period 5 and giving rise to 

output functions of period 5 itself. Obviously, only coefficients 

r k 7 3 6 9 

corresponding to eigensequences 4 e j for e = 1, ct , a , a and 
a 12 appear at the input and at output, since they are the 
sequences of period 5. In each case it is seen that 
(i) Y * S(e) U 


for all e 
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Table 4.5 


Spectral Coefficients for Input and Output for 
Sequences of Period 5; 2 -St age Scrambler 


t 


Critical 

States X 
c 

X 1 X 2 


Input/ 

Output 

Sequences 


x 


Input/Output Spectral Coefficients* 



10000,... 

1 

0 

0 

1 

0 

1 

0 

0 

1 

0 

1 

0 

0 

1 0 


01101,... 

1 

0 

0 

14 

0 

12 

0 

0 

15 

0 

8 

0 

0 


01000,... 

1 

0 

0 

13 

0 

10 

0 

0 

7 

0 

4 

0 

0 

o . 1 


10110,... 

1 

0 

0 

11 

0 

6 

0 

0 

6 

0 

ii 

0 

0 


00100,... 

1 

0 

0 

10 

0 

4 

0 

0 

13 

0 

7 

0 

0 

1 1 

" 

01011,... 

1 

0 

0 

8 

0 

15 

0 

0 

12 

0 

14 

0 

0 


00010, . . . 

1 

0 

0 

7 

0 

13 

0 

0 

4 

0 

10 

0 

0 

1 0 


10101,... 

1 

0 

0 

5 

0 

9 

0 

0 

3 

0 

2 

0 

0 


00001,... 

1 

0 

0 

4 

0 

7 

0 

0 

10 

0 

13 

0 

0 

0 1 


11010,... 

1 

0 

0 

2 

0 

3 

0 

0 

9 

0 

5 

0 

0 

* Mapping of 

elements: 

0 — 

0 

/ 

i 

* : 

1 

* a 1 , 

i 

*5* | 

3/ 1 

/ • • 

.,14 
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(ii) U X(e) 

e 1 e 


+ X 


* 0 


where X^ is the critical state for the sequence and e = 1 , cc ', 

6 9 12 

a * cc < and a . This has to be so, otherwise output period will 
be 1CM (15, 3) = 15. >* 


From Section 3.9 it follows that the above conditions are satis- 
fied in general for any sequence of period five giving rise to 
output sequence of period five. 

Now let us consider the case when the input period of 15 
resulting in output period of 15, irrespective of initial state. 
The cases of input period 5/3 resulting in output period 15 is 
also covered in this. Spectral coefficients at the input and at 
output for 13 basis function. When the initial state of the 
scrambler is zero is given in Table 4.6. The sequence 
111000000000000,111... and its shifted versions are taken as 


basis functions because. 


Let the input transform be 


■P-ill- 
,15 


1 + d 

For the two stage scrambler from Section 3.5, the transfer 
function is given by 


H(d) 


d + d + 1 

Therefore the output transform is given by 


¥<d) 


mi 


Tc * 9 * 

1 + d (1 + d + d ) 


deg p(d) < 14 


In the denominator of Y(d> there are two factors of (1 + d + d ) 
Hence the exponent of the denominator can be 15 only if p(d) 



has a factor (1 + d + d 2 ). in that case only the output period 
will be 15 irrespective of initial state. Otherwise the output 
period will be 30. Hence p(d) can be anything of the kind given 

t>y 

p(d) = (1 + d + d 2 ) (a Q + ajd 1 + a 2 d 2 + ... + a ]c d k ) , k 4 12 


Hence it follows that any function resulting in an output period 
of 15 should be some linear combination of these 13 sequences 
given in the table. Last three rows gives the zero input 
response for different states and their spectral contents. 

It can be seen that 


(i) 

(ii) 

(iii) 


= S(e) U , e £• (e)| eigenvalue is defined for 

The new coefficients are due to the 2 ero input response 
of the state X T given by 


*T 



X(e) 


for all e. 


The spectral coefficients at the input corresponding to 


those places where eigenvalue is not defined are zero. 

In general any output sequence is some linear combination 
of these 13 output sequences, plus any one of the last three 
output sequences for non-zero initial state. Hence it follows 
that for any such sequence the conditions (i) , (ii) and (iii) 
are satisfied and the state Xp for any non— zero initial state 

X is given by 

ct 
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(d.) Eigensequences from GF(2 )j in this extension field, we 
will consider cases having input and output periods 31. Table 
4.7 gives the eigenvalues and unique state for each eigensequences. 
Spectral coefficients for 31 basis functions at input and at 
output along with their critical states are given in Table 4.8. 

It can be seen that 


(i) Y e = U e S(e) ' 


(ii) X(e) = 

q e c 


for all e. 


in each case. Hence from Section 3.9 it follows that in all 
cases of input period of 31 resulting in an output period of 31, 
the conditions (i) and (ii) are satisfied. 


4.7.2 Three Stage Scrambler: 

We will analyse the three stage scrambler structure 
shown in Fig. 4.5. for the following cases 



Figure 4.5 Three-stage Scrambler; c(x)-x +x + l 
(i) Input and output periods 3 
(ii) Input and output period 7 
(iii) Input and output period 5 
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(iv) Input and output period 15 

(v) Input and output period 31 

(vi) Input period 3/5 , output period 15. 

2 

(a) Extension field GF(2 ): The eigenvalues and corresponding 

unique states are listed in Table 4,9, 


Table 4.9 Eigenvalues and Unique States in GF(2 2 ); 

3 • 

e(x)=x +x+l 


e 

p. 

1 

1 

« 

$ 

I 

1 

1 

1 

X 1 

X 

e 

X 2 

1 

( 

1 

1 

1 

X 3 I 

s(e> 

i 


1 

1 

X 

1 

a 


a 

1 

a 2 

„2 
c t 

« 2 


a 2 

1 

a 

a 


The spectral coefficients at input and output for three basis 
functions , given as Input when the scrambler is in critical 
state is listed in Table 4.10. From the table it can be seen 
that 

(i) Y e - S(e) U e 

(11) £> e X(e) = X c 

e 

From Section 3.9, it follows that in all cases where input and 

i 

output periods are three, above two properties hold good. Note 

, . 4 

that the analysis for period 3 can be done in GF(2 ) also since 

GF(2 2 ) C GF(2 4 ). 
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Table 4*10 Spectral Coefficients at Input and Output for 
Sequences of Period 3; c(x) >= x 3 + x + 1 


Critical 

State X 
c 

X 1 *2 X 3 

Input/Output Sequence 

Input/Output Spectral 
Coefficients 



100-,.*. 

1 

1 

1 

0 1 

0 



o 




010,... 

1 

a 2 

a 


010 , . . . 

1 

a 2 

a 

1 0 

0 




0 



001,.. . 

1 

a 

a 



001 , . . . 

1 

a 

2 

a 

0 0 

1 







100, ... 

1 

1 

1 


3 

(b) Eigensequences from GF(2 ): 

3 

Eigenvalues for eigenfunctions in GF(2 ) and corresponding 


unique initial states are 

given in Table 4.11 . 


Table 4.11 

Eigenvalues 
c(x) = x 3 + 

and Unique 
x + 1 

States in 

GF(2 3 ) 

e 

1 

l 

1 

x e 

t 1 

1 

1 

1 

s(e). 


| X 1 

J ....... 

X 2 

X 3 j 

- - 1 - 


i 

1 

1 

i 

1 

3 

a 

2 

a 

a 6 

a 3 

a 5 

5 

a 

a 

a 3 

a 5 

.a & 

6 

a 

4 

a 

5 

a 

a 6 

a 3 


2 4 

Notes Eigenvalues not defined for e = cc, a and a . 
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Prom Equation (4.23) we have basis function for the 
sequence space# any Sequence from which results in output period 
of 7 irrespective of initial state. Any other sequence of period 
7 outside this space gives rise to output period 14. Spectral 
coefficients for these basis functions at input and at output 
when the initial state is zero is listed in Table 4.12. In the 
last row of the table is given# spectral coefficients for the 
zero input response when initial state is 001. Since for other 
initial states# zero input response is some shifted version of 
the same# spectral coefficients exist exactly at the same loca- 
tions, but their coefficients are different. 


Table 4.12 Spectral Coefficients at Input and Output for 
Sequences of Period 7; c(x) = x 3 + x + 1 



Initial 


« 

J Input/Output 

• 

1 

J Input/Output 

Spectral Coefficients 


State X„ 

J Sequence 

; for 

e - 






X 1 

c 

*2 X 3 

• 

1 

1 

1 

• * 

1 

i 1 

1 




2 

a 

3 

a 

4 

a 

5 

a 

6 

a 




1111111#... 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1101000#. .. 

1 

a 

a 2 

0 

a 4 

0 

0 




1110100#... 

0 

0 

0 

1 

0 

1 

1 

0 

0 

0 

1100000# . . . 

0 

a 2 

a 4 

a 5 

a 


a 3 




0111010# . . . 

0 

0 

0 

a 4 

0 

v 

a 

0 

0 

0 

0110000#. . . 

0 

a 

a 2 

a 2 

4 

a 

a 

4 

a 




0011101#... 

0 

0 

0 

a 

0 

; 4 ~ 

i 

p . 
M 

1 

0 

0 

0 

0011000#... 

0 

1 

1 

J> 

a 

1 

a 3 

5 

a 




0000000# . . . 

0 

0 

0 

i 

O 1 

1 

0 

o 

1 

0 

0 

0 

1 

1011100#... 

0 

a 3 

a 6 

0 

£ 

0 

0 
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Prom Table 4.12 and linearity o£ the scrambler it follows that 
for any sequence of period 7 and giving rise to output period 7, 


(i) 

(ii) 

(ill) 


Y Q « U e S(e) / e q (e j eigenvalue is defined for| 
Spectral coefficients at locations where eigenvalues are 
not defined are zero at the input side. 

U X(e) + X_ generates sequences which give rise to 
new spectral coefficients at the output , at locations 
where eigenvalues are not defined. 


(c) Eigenfunctions from GF(2 4 ): Eigenvalues and corresponding 

4 

unique states for all eigenfunctions in GF(2 ) are. listed in 
Table 4.13. Let us first consider the case of input and output 
sequence both having periodicity 5. In Table 4.14 we have 
listed the spectral coefficients for five basis functions at 
input and at output when the scrambler is in critical state. In 
all cases, it can be seen that 
(i) Y a = u s(e) , for all e 

(ii) u e X<e> = x c . % 

e 

This has to be so otherwise the period will be LCM (5, 7), 
i.e. 35. 


From Section 3.9 it follows that the above two properties are 
satisfied in all cases where the input sequence of period 5 
results in output sequence of period 5 itself. 

New let us consider an input sequence of period 15 and 
giving rise to output period of 15 itself. Table 4.15 gives 
the spectral coefficients of 15 basis input sequences and 
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Table 4,13 

Eigenvalues 
c(x) « x3 + 

and Unique 
x + 1 

States 

in GF(2 4 ); 


e 

f r r™"*™*' ' r ' n rTr 1 1 

i 

» 

X(e) 


• 

! S(e) 



|: x i 

i 

f . . . 

x 2 

X 3 



i 

i 

1 

1 

1 


a 

a 10 

a 9 

a 8 

a 11 


a 2 

a 5 

a 3 

a 

a 7 


a 3 

a 2 

a^ 4 


a 5 


a 4 

a 10 

c£ 

a 2 

a 44 


a 5 

a 8 

1 

dLO 

a 10' 


a 6 

d 4 

a 13 

a 7 



a 7 

ct 9 

a? 

a 10 

a 


« 8 

« 5 

a 12 

a 4 

a 13 


a 9 

a 

a 7 

a 3 - 3 

a 


a 10 

a 10 

l" 

a 3 

a 5 


d 11 

a 12 

a 

■ 

a 8 


a 12 

a 8 

o ^ 

a 14 

a 5 


13 

a 

a 6 

cP 

a 10 

a 4 


a 14 

a 3 

$ 

a 8 

a 2 



of the resulting output sequences when the scrambler is in 
critical state. It can be seen that in each case 
(i) « S(e) U , for all e 

(ii> 21! u e x < e > - * c • 

0 ’ . 

This has to happen otherwise period will be l/3M(15.,7), i.e. 105. 
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Table 4,14 Spectral Coefficients at Input and Output for 
Sequence of Period 5; c(x) = x 3 + x + 1 


Critical 

states X 
c 


Input/ 

Output 

Sequences 


Input/Output Spectral Coefficients* 



» 

10000, . . . 

1 

0 

0 

1 

0 

0 

1 

0 

0 

1 

0 

0 

1 

0 

0 

100 


11001,... 

1 

0 

0 

6 

0 

0 

11 

0 

0 

11 

0 

0 

6 

0 

0 

«*. — 

» mm mm .mm ^ mm mm 

- 

— 

— 

— 

— 

— - 

— 

— — 

— 

— 

- 

- - 

- - 

— 

— 

001 

01000,... 

1 

0 

0 

13 

0 

0 

10 

0 

0 

7 

0 

0 

4 

0 

0 


11100,.., 

1 

0 

0 

3 

0 

0 

5 

0 

0 

3 

0 

0 

9 

0 

0 

mm mm mm 

— 

- 

— 

— 

— 

— 

— 

— 

— 

— 

— 

- 

— 

— 

— 


00100,.,. 

1 

0 

0 

10 

0 

0 

4 

0 

0 

13 

0 

0 

7 

0 

0 

Oil 


OHIO,... 

1 

0 

0 

15 

0 

0 

14 

0 

0 

8 

0 

0 

12 

0 

0 

1 1 1 

00010, . . . 

1 

0 

0 

7 

0 

0 

13 

0 

0 

4 

0 

0 

10 

0 

0 

X x x 

00111,... 

1 

0 

0 

12 

0 

0 

8 

0 

0 

14 

0 

0 

15 

0 

0 


00001 , . . . 

1 

0 

0 

4 

0 

0 

7 

0 

0 

10 

0 

0 

13 

0 

0 

110 


10011,... 

1 

0 

0 

9 

0 

0 

2 

0 

0 

5 

0 

0 

3 

0 

0 


*Mapping of elements as in Table 4.5 
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From Section 3.9, it follows that in all cases where an input 
sequence of period 15 giving rise to output period 15, the 
above two conditions are satisfied. The case of input period 
3/5 and output period 3/5 is also covered in this. 

(b) Eigensequences from GF(2 ): Eigenvalues ' and unique states 

for all eigensequences are given in Table 4. Id., We will consider 

the case of input and output sequences having period 31. The 

spectral coefficients for 31 basis functions at input and for 

the resulting output sequences when the scrambler is in critical 

state are given in Table 4.17. It can be seen that 

(i) Y « U S(e) , for all e 
6 © 

<ii) X" u e X(e) = X c 

e .. 

in each case. From Section 3.9 it follows that any sequence of 
input period 31 resulting in a sequence of output period 31 
satisfies the conditions (i) and (ii) « 

Exactly similar properties are observed in the case of 

the three stage scrambler having characteristic polynomial 
3 2 

C(x) =? x + x +1. This can be verified from Tables B.l to B.6 
(Appendix B) . Since the cases of input periods 3 and 5 resulting 
in output periods of 3 and 5 respectively are covered in the 
case of input and output period 15, separate tables for those 
cases are not attached. 
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The spectral coefficients at input and output for three basis 
functions , given as input when the scrambler is in critical 
state is listed in Table 4.10. Prom the table it can be seen 

that 
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From Equation (4.23) we have basis function for the 
sequence space,, any sequence from which results in output period 
^ irrespective of initial state. Any other sequence of period 
7 outside this space gives rise to output period 14. Spectral 
coefficients for these basis functions at input and at output 
when the initial state is zero is listed in Table 4.12. In the 
last row of the table is given, spectral coefficients for the 
zero input response when initial state is 001. Since for other 
initial states, zero input response is some shifted version of 
the same, spectral coefficients exist exactly at the same loca- 
tions, but their coefficients are different. 


Table 4.12 Spectral Coefficients at Input and Output for 
Sequences of Period 7? c(x) = x 3 + x + 1 


Initial 

State x c 

*1 *2 X 3 

1 

Input/Out put 
Sequence 

5 Input /Out put Spectral Coefficients 
I for e = 

* 2 3 4 S 6 

.1 a a a a a 




1111111,... 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1101000,... 

1 

a 

a 2 

i 

1 o 

1 

a 4 

I 

1 o 

I 

0 




1110100,... 

0 

0 

0 

1 

0 

1 

1 

0 

0 

0 

1100000,... 

0 

a 2 

a 4 

a 5 

a 

& 

c? 



m- mum mm . , 

0111010,.. . 

0 

0 

0 

>' 

0 


a 

0 

0 

0 

0110000,... 

0 

a 

a 2 

a 2 

A 

a 

a 

a 4 

I 

1 

mm mmm m 


0011101,... 

0 

0 

0 

a 

0 

■; 4 ~ 

V' 

0 

0 

0 

0011000, . . . 

0 

1 

1 

6 

a 

1 

0? 

a 5 


«mm «M» m 

jjjBpfeil 

0000000, . . . 

0 

0 

0 

0 

0 

o 

'i 

0 

0 

0 


1011100,... 

0 

a 3 

a 6 

0 

2 

0 

0 
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From Table 4.12 and linearity of the scrambler it follows that 
for any sequence -of period 7 and giving rise to output period 7 , 

^ Y e ~ u e S(e) , e £ (e j eigenvalue is defined for £e k |) 
Spectral coefficients at locations where eigenvalues are 
not defined are zero at the input side. 

-<—-4 u e (ii) * * * * * * * x (e) + X & generates sequences which give rise to 
new spectral coefficients at the output, at locations 
where eigenvalues are not defined. 

(c) Eigenfunctions from GF(2 4 ): Eigenvalues and corresponding 

unique states for all eigenfunctions in GF(2 4 ) are listed in 
Table 4.13. Let us first consider the case of input and output 
sequence both having periodicity 5. in Table 4.14 we have 
listed the spectral coefficients for five basis functions at 
input and at output when the scrambler is in critical state. In 
all cases , it can be seen that 
(i) Y Q = U e S(e) , for all e 

(ii) Z* U X(e) = X, 
e c 

This has to be so otherwise the period will be LCM (5, 7), 

i.e. 35. 

From Section 3.9 it follows that the above two properties are 

satisfied in all cases where the input sequence of period 5 

results in output sequence of period 5 itself. 

Now let us consider an input sequence of period 15 and 

giving rise to output period of 15 itself. Table 4.15 gives 

the spectral coefficients of 15 basis input sequences and 
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© 4.14 Spectral Coefficients at Input and Output for 
Sequence of Period 5; c(x) = x 3 + x + 1 


Critical 
States X 



Input/Output Spectral Coefficients* 


100 

10000 1 

0 

0 

1 

0 

0 

1 

0 

0 

1 

0 

0 

1 

0 

0 


11001,... 1 

0 

0 

6 

0 

0 

11 

0 

0 

11 

0 

0 

6 

0 

0 

001 

01000,... 1 

0 

0 

13 

0 

0 

1 O 

0 

0 

7 

0 

0 

4 

0 

0 

11100,... 1 

0 

0 

3 

0 

0 

5 

0 

0 

3 

0 

0 

9 

0 

0 

oil 

00100,... 1 

0 

0 

10 

0 

0 

4 

0 

0 

13 

0 

0 

7 

0 

0 

OHIO,... 1 

0 

0 

15 

0 

0 

14 

0 

0 

8 

0 

0 

12 

0 

0 

111 

00010,... 1 

0 

0 

7 

0 

0 

13 

•» 

0 

0 

4 

0 

0 

10 

0 

0 

00111,... 1 

0 

0 

12 

0 

0 

8 

0 

0 

14 

0 

0 

15 

0 

0 

110 

00001,... 1 

0 

0 

4 

0 

0 

7 

0 

0 

10 

0 

0 

13 

0 

0 


10011 ,... 1 

0 

0 

9 

0 

0 

2 

0 

0 

5 

0 

0 

3 

0 

0 


*Mapping of elements as in Table 4.5. 






with following mapping for input/output spectral coefficients column 
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From Section 3 q i *. ^ , , 

' follcws that In all cases where an input 

SeqUenCS ° f Perl ° d 15 ^ng rise to output period IS, the 
above two conditions are satisfied. The case of input period 

3/S and output period 3/5 is also covered in this. 

<b> Eigensequences from GFC 2 5 ), Eigenvalues and unique states 
for all eigensequences are given in Table 4 . 16 .. we will consider 
se of input and output sequences having period 31. The 
spectral coefficients for 31 basis functions at input and for 
the resulting output sequences when the scrambler is in critical 
state are given in Table 4.17. It can be seen that 
<i) Y e = u e s(e> , for all e 


(li > U e X(e) = 

e .. 


X 


In each case. From Section 3.9 it follows that any sequence of 
input period 31 resulting in a sequence of output period 31 
satisfies the conditions (i) and (ii). 


Exactly similar properties are observed in the case of 

the three stage scrambler having characteristic polynomial 

3 2 

C(x) = x + x +1. This can be verified from Tables B.l to B.6 
(Appendix B) . Since the cases of input periods 3 and 5 resulting 
in output periods of 3 and 5 respectively are covered in the 
case of input and output period 15, separate tables for those 
cases are not attached. 
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*7*3 Four Stage Scramblers: 

characteristic « “* T ^ *- 

tic polynomial C(x) = x 4 4 . v . . 

+ x+ 1 • shown in Figure 4.6 



Input 


Output 

4 ’ 6 Scrambler; c<x> 

^ tWS 0386 We W1U do s Pcctral analysis which covers 
the following cases. 

U) input period 3 output period 3 
<ii) Input period 7 output period 7 
(ill) Input period 5 output period 5 
Uv) In Put period 3/5 output period 15 
(v) Input period 15 output period 15 
(vi) Input period 31 output period 31. 

(a) Exponent function from GF(2 2 ): Eigenvalues and unique 
States for all eigensequences are given in Table 4.18. Table 
4.19 gives input and output spectral coefficients for 3 basis 
sequences of period 3 and corresponding critical states. p rom 
Section 3.9 and these tables it follows that 


** ^ * le ■ . 


:(e) = x 





# 
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for any sequence of input and output periodicity 3. This is 
justified by Equation ( 4 . 50 ). 


ble 4.18 ^envalues and Unique States in GF(2 2 ); 
c ' x ' - x* + x + 1 



Table 4.19 Spectral Coefficients at Input and Output for 
Sequences of Period 3; c(x) = x 4 + x +1 


Critical State X c J Input/Output j Input/Output Spectral 
^ X 2 . x 3 x 4 ; Sequences ; Coefficients 


0 0 


100 , . . 

001 ,.. 


1 

1 


1 

a 


a 


0 


0 


010 ,... 

100 ,*;. 


001 ,... 

010 ,... 


— — I 


1 

1 


1 

1 


a 

1 


a 


a 


a 

1 

a 2 

a 
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(b) Eigensequences f r 
for all elgensequences 
spectral coefficients ■ 
critical states are gi- 
from the table it can ] 


.20. Input and output 
of period 7 with their 


ler the case of input and output periods 
ifficients at input and output for 5 basis 
along with their critical states are 
From Section 3.9 and from this table it 


2 we have obtained a set of basis functions which 
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consisting of such sequences. The spectral coefficients at 
input and output for these basis functions when the scrambler 
is in zero state are listed in Table 4.24. The last row in the 
table gives the spectral coefficients of zero input response 
when the scrambler state is 0001. The zero input responses for 
other states have spectral contents similar, only coefficients 
being different. In each case it can be seen that 
<i) Y q * U e S(e) , e ( e | eigenvalue is defined for (e^J) 

(ii) The spectral coefficients at input in those places where 
eigenvalue is not defined are zero. 

(iii) The new spectral coefficients appearing at output are 

due to the zero input response for the state ^ given by 

*T = 2>e x(e > + x a 



0 4 ' 21 values and Unique states in GP( 2 3 ),- 


Critical State X 


c ; 

X • Sequences 


Input/Output Spectral 



1000000,... 

1 

1 

1 

1 

1 

1 

i 

1011110,... 

1 


a 3 

£ 


a 3 

a 5 

0100000,... 

1 


a 5 

1 

i 

y 

i 

R to 1 

a 

0101111,... 

1 

a 4 

a 

a 3 


a 5 

3 

0010000,... 

1 

a 5 


a 

~y 

1 

| 

a 2 

1010111,... 

1 

a 3 

a 6 

i 

3 

1 

1 

0001000,... 

1 

a 4 

a 

<3 

c? 

a 6 

Q 3 

1101011,... 

1 

a 2 

4 

a 

a 4 

a 

3 

a 

0000100,... 

1 

« 3 ” 

► mm mm:. 

a 6 

y 

V 

a 

a 4 

1110101 

1 

a 

a 2 

a 

a 4 

a 4 

a 2 

0000010,... 

1 

a 2 

a 4 

<£ 

a 

3 ' 

1 

1111010,... 

1 

1 

1 

5 

a 

1 

a 6 

c? 

0000001,... 

1 

a 

y 

a 3 

" .mm,: ■■■» 

£ 

7 a 6 

0111101,... 

1 

a 6 

a 5 

u 2 

£ 

a 

4 

a 
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Table 4.22 


Eigenvalues and ~ 
c(x) = x 4 + x + 1 


States in 


) ; 


e 

t 

1 

1 

■, f 


X(e) 


f 

• 

t 

t 

S(e) 



i x i 

— 1 

X 2 

X 3 

X 4 

l 

f 

* 

I 



1 

1 

1 

1 

1 


1 


a 3 

a 4 

a 

a 13 

a 10 


a 7 



1 


a? 

1 


a 5 


a 6 

a8 

a 2 


a 5 


a 14 


a 7 

a 11 

& 

a 1 - 2 

a 5 


a 3 


a 9 

a 2 

8 

a 

14 

a 

a 5 


11 

a 


a 10 

1 

a 5 

10 

a 

1 


10 

a 


a 11 

a 13 


£ 

a 10 


0 ? 


12 

a 

a 

4 

a 

a 7 

10 

a 


13 

a 


13 

a 

14 

a 

a 

3 

a 

5 

a 


12 

a 


14 

a 

7 

a 

8 

a 

9 

a 

10 

a 


6 

a 



where X is the actual state and e < 0 . (e | eigenvalue is 
■ k 

defined for £e j ) . 

Prom the linearity of the scrambler it follows that the entire 
space of sequences having input and output periods 15 irrespec- 
tive of initial state satisfies the conditions (i), (ii) and 

(iii) . 
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Input/ 

Output 

Sequences 


10000 


1101 


01011 


01000 


1010 


10101 


00100 


0101 


11010 


00010 


01101 


00001 


0110 


10110 


Mapping as in Table 4.5 


The eigenvalues with their 


corresponding unique states are given in Table B.7. The spectral 
coefficients and output for 31 basis functions along with 

their critical states are listed in Table B.8. From this table it 


From Section 3.9 it follows that these two conditions are 
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4.7.4 Five Stage Scramblers: 

The five stage scrambler having the characteristic 
polynomial c(x) = x +x + 1 is given in Figure 4.7. In this 



Output 


Figure 4.7 Five-stage Scrambler; c(x) = x^ 4 - x^ +1. 

we will analyse the scrambler for the f ol lowing cases. 

(i) Input period 7 output period 7 

(ii) Input periods 3, 5 and 15 giving rise to output periods 
3, 5 and 15 respectively 

(iii) Input period 31 and output period 31. 

3 

(a) Eigensequences from GF(2 ): The eigenvalues and their 

corresponding unique states are listed in Table 4.25. The 
spectral coefficients at input and at output for 7 basis func- 
tions of period 7 and their critical states are listed in Table 
4 . 26,. From this table and from Section 3.9 it follows that for 
any sequence of input and output periodicity 7, 

(i) Y = U S(e) / for all e 
e e 

Ui> u x(e)' = x o . 

e 

This is justified by Equation (4.50). 




1 

1 

1 

1 

1 

. 1 

1 


3 

2 



6 

4 

a 

a 

a 

a 

1 

a 

a 

€? 
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a 4 

a? 

1 

a 5 

a 

a 3 

a 

a 5 

a 2 - 

a 6 

a 3 

<# 

a 4 

a 5 

a 

a 4 

1 

0? 

€? 

a 5 

0^ 

a 6 

a 

a 3 

a 5 

C? 

a 6 

a 2 

a 3 

a 4 

a 5 


a 






Table 4.26 Spectral Coefficients at Input and Output for 
Sequence of Period 7; c(x) = x 5 + x 2 + 1 


Critical | 

States X [ Sequences 


Input/Output Spectral 
Coefficients 


00111 

1000000,... . 

* 

1 

1 

i 

1 

1 

1 

i 


1111100, ... 

1 

„4 

a 

a 

a 4 

a? 

£ 

a 


0100000, .. . 

1 


’a 5 ” 

1 

1 

1 

y 

i 

1 (Nfcj 
1 

a 

0111 


a 6 

a 4 



0111110,... 

1 

a 3 

a 


a 2 

11111 

0010000,... 

1 

5 

a 

a 3 

a 

6 

a 

£ 

a 2 

a 3 

0011111, ... 

1 

a 2 

a 4 

«? 

a 

a 6 

11110 

0001000, . . . 

1 

"a 4 " 

a 

~y 

a 2 

a 6 

a 3 

1001111, ... 

1 

a 

a 2 

£ 

a 4 

a 

a 4 


0000100,. .. 

1 

3 

a 

6 

a 

~y 

’y 

a 

a 4 

11100 

1100111, ... 

1 

1 

1 

a 6 

i 




0000010, . . . 

1 


a 4 

6 

a 

a 

£ 

~£ 

11001 

1110011,... 


6 

a 

5 

3 

3 

5 

6 


1 

a 

a 

or 

a 

a 

1 

1 

1 

1 

1 

1 



• — 

- — — 

’ — - 

- — 


- — — 


0000001,... 

1 

a 

a 2 

a 3 

a 4 

5 

a 

6 

a 


10011 


1111001 ,.., 


a 


a 


6 


a 
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(b) Eigensequences from GF(2 4 ): The eigenvalues and corres- 

ponding unique states are listed in Table 4.27. The spectral 
coefficients at input and at output for 15 basis functions of 
period 15 and their critical states are listed in Table 4.28. 
From this table and from Section 3.9 it follows that for any 
sequence of input and output periodicity 15 , 

(1) ' Y e = u e S(e) ' for all e 


(ii) T* U x(e> = X 

-^-t e c 


This is j list i fled by Equation (4.50). 

This covers the case of input and output periodicities 3 and 5 
also. 


(c) Eigensequences from GF ( 2 ): Eigenvalues and corresponding 

unique states are listed in Table 4.29. The eigenvalues are not 
defined for e = a, a , a 4 / a® and Let us now consider the 

case of input sequences of period 31 giving rise to output 
sequences of period 31 irrespective of initial state of the 
scrambler. Preceeding as in Illustration 4.1, a set of 26 basis 
functions which span the vector space consisting- of such sequ- 
ences are obtained. The spectral coefficients at input and at 
output for these basis functions when the scrambler is in zero 
state are listed in Table 4.30. The last row in the table gives 
the spectral coefficient of zero input response when the scram- 
bler state is 00001. The zero input responses for other states 
have spectral contents similar, only coefficients are different. 
In each case it can be seen that 



153 


Table 4.27* 

Eigenvalues and Unique 
c(x}' * x^ + x^ +1 

* 

States 

in GF(2 4 ); 


e 


*2 

X<e) 

X 3 

X„ 

4 

Xsl 

1 

S(e); 

i 

1 

1 

1 

1 

1 

1 

2 

1 

15 

14 

13 

12 

2 

3 

1 

14 

12 

10. 

8 

3 

4 

7 

4 

1 

13 

10 

1G 

5 

1 

12 

8 

4 

15 

5 

6 

6 

1 

11 

6 

1 

11 


7 13 7 1 IQ 4 4 

8 8 11 4 12 5 IQ 


:l*lg|f j 


9 

1 

8 

15 

7 

sigiiliiig 

14 

9 

10 

4 

10 

1 

7 

13 

13 

11 

11 

1 

6 

11 

1 

6 

12 

2 

6 

10 

14 

3 

13 

13 

10 

13 

1 

4 

7 

7 

14 , 

9 

11 

13 

15 

2 

7 

15 

5 

6 

7 

8 

9 

4 
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(iii)' The new 


where X a is the initial state of the scrambler. 

Prom the linearity of the scrambler it follows that the entire 
a P a< - e sequences having input and. output periods equal to 31 

irrespective of initial state satisfies the conditions (i)‘, (ii > 
and (iii):. 


5-stage scramblers. It can be verified from the corresponding 
Tables B.l8 to B.23 for the scrambler with c(x) = x^ 


From the analysis carried out it is seen that 
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Y e ~ U e s k%<= (e ( eigenvalue is defined f or £ e X j ), 
Spectral coefficients at input, of those eigensequences 
for which eigenvalues are not defined, are zero, and 
The new coefficients which appear at the output are 
due to the zero input response of the scrairibler for 
the state ^ given by 

*T = Z! U e X(e); + x a 

€3 


where X a is the actual initial state of the scrambler. 



CHAPTER 5 


CONCLUSION 

In .this chapter we summarize various results obtained 
and observations made in this thesis and make suggestions for 
further work in this area. 

5.1. Conclusion 

The important results obtained concerning the periodic 

output of the scrambler when excited by a periodic input are: 

r ■■ ■ ■ : • • ' i 1 ), " , r > mswmSm if 

1. If the periods of input sequence is not a multiple of the 
cycle length q of the scrambler, then but for a unique 
initial state of the scrambler for which the period remains 
unaltered, the period of the output sequence is LCM(s, q). . 

2. If the period s of the input sequence is a multiple of the 
cycle length of the scrambler, say s is seme rq, then the 
period of the output sequence is also a multiple of the 
cycle length i.e. the period is some r*q (r* not necess- 
arily different from r);. 

3. There is no advantage in scrambling a M-sequence of length 
equal to cycle length of the scrambler and generated by a 
structure different from that of the scrambler. 

From the spectral analysis of scramblers, the important 
observations made are: 

For a n-stage scrambler eigenvalues for n eigensequences 
over an extension field GF(2 n ) are not defined. However, 


1 


161 


eigenvalues are defined « n *. ,n , nu 

erinea tor all eigensequences over GF(2 > 

when n is not a divisor of m. 

2. For input and output sequences of period not equal to cycle 
length of the scrambler, 

(a) Y q = u e s(eX for all e 

and (b > X! u e x(e > = X 

e c 

3. For input and output sequences of period, a multiple of 
cycle length, 

(a) Y e * u e S(eX ,e(J (e | eigenvalue is defined for ^e*}) 

(b) spectral coefficients at input, of those eigensequences 
for which eigenvalues are not defined, are . * z er o, 

and 

(c ) ' the new spectral coefficients which appear at the 

output are due to the zero input response of the 
scrambler for the state X^, given by 



where X a is the initial state of the scrambler. 

5.2. Suggesticns , for Further Research 

Performance of scramblers defined over any modular field 
GF ( p )' in general can be studied using techniques similar to those 
employed in Chapter 3. 

In this thesis spectral theory for LSCs defined over a 
modular field is used to carry out the spectral analysis of the 
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input and output sequences and observations as outlined in 
Section 5.1 are made based on this analysis. Hie analysis 
reported in the thesis holds good for input and output sequences 
with periodicities equal to 2 n -l or to a divisor of 2 n -l for 
some integer n. Further investigation may be carried out to see 
whether LSCs can be studied completely in a manner analogous to 
frequency domain analysis of continuous linear time-invariant 
systems. In other words whether it is possible to develop a 
theory similar to that of Fourier transform, by exploiting 

the concept of formal derivitive discussed in Chapter 3. 
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APPENDIX A 


OPERATION TABLES for GW2 2 ), GF(2 3 ). GP(2 4 ) AND GF(2 5 


1 


This appendix presents the operation tables for the 
Galois fields GF(2 ), for n = 2, 3, 4 and 5. The irreducible 
polynomial used to generate these fields is also indicated. 

1. GF< 2 2 ); ( 0 , i, a, a 2 ): g< x ) = x 2 + x + 1 

Table A.l Addition Table for GP(2 2 ) 


0 

1 

a 


0 

0 

1 

a 


1 

1 

0 

d 

a 

a 


a 

0 

1 


a 

a 

1 

0 


Multiplication: a* . a^ = a^ i+3 * I mod 3 


2. GF(2 3 )? (0, 1/ a, a 2 ,...,a^); g(x) = x 3 + x+ 1 



Table A, 2 Addition Table for GF(2 3 ) 


+ 

0 

1 

a 

a 2 

a 3 

a 4 

5 

a 

a 6 

0 

0 

1 

a 

a 2 

a 3 

4 

a 

5 

a 

6 

a 

1 

1 

0 

3 

a 

a 6 

a 

a 5 

4 

a 

a 2 

a 

a 

a 3 

0 

4 

a 

1 

a 2 

6 

a 

S 

a 

a? 

a 2 

a 6 

a 4 

0 

5 

a 

a 

3 

a 

1 

3 

a 

3 

a 

a 

1 

5 

a 

0 

6 

a 

a 2 

a 4 

4 

a 

a 4 

a 5 

2 

a 

a 

6 

a 

0 

1 

a 3 

5 

a 

a 5 

4 

a 

6 

a 

3 

a 

a 2 

1 

0 

0 ~ 

6 

a 

6 

a 

a 2 

5 

a 

1 

4 

a 

3 

a 

a 

0 

Multiplication : 

a 1 . 

a-* sr 

(i+j ) J mod 7 
a 





9H 



(0, 1, a. 


2 <| ^ 

a *•••*&■ )/ g(x) = a 4 + x + 1 


3. GP(2 4 ); 

Table A. 3* Addition Table for GF(2 4 ) 
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1 
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7 

8 
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10 

11 

12 

13 

14 

15 
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11 
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13 

14 

15 
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15 
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10 
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13 
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7 

4 
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6 

10 
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15 

11 
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13 
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3 
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11 
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13 
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12 
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10 
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15 

14 
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4 

15 

10 
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12 
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14 
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13 
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11 
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5 
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11 
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9 

13 

4 
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15 
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14 
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12 

10 

6 
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11 
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2 

12 
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10 

14 
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15 
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13 
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14 

12 
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13 
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15 
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14 
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14 
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10 
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id 

10 
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